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Potential  Advantages  of  the  Spaced  Antenna  Method  for 
Operational  Wind  Profiling 
J.  Rottger1  and  M.  P.  Larsen2 
1 .  Introduction 

The  problem  of  very  short-range  forecasting  is  twofold.  It  is  necessary  to 
understand  the  processes  that  are  being  forecasted,  and  data  appropriate  to  the  scale  of 
interest  has  to  be  generated  .  Coherent  VHF  and  tJHF  radars  are  being  used  for 
operational  wind  profiling  and  are  providing  part  of  the  solution  to  the  data 
acquisition  problem  (Larsen,  1983)*  The  Profiler  system  operated  by  the  Wave 
Propagation  Laboratory  at  MOAA  has  already  shown  great  promise  (Strauch,  1981;  Strauch 
et  al. ,  1982).  As  a  result,  plans  are  being  considered  for  expanding  the  network  of 
radars  to  cover  a  larger  area  of  the  country. 

The  Profiler  uses  what  is  commonly  referred  to  as  the  Doppler  method  for 
measuring  winds.  Two  beams  are  pointed  off-vertical,  and  the  Doppler  shift  of  the  echo 
determines  the  llne-of-sight  velocity.  The  velocity  components  along  the  beams  are  then 
translated  to  horizontal  wind  components.  While  there  is  no  doubt  that  the  Doppler 
method  is  adequate  for  wind  profiling,  we  want  to  discuss  a  number  of  possible 
advantages  of  the  spaced  antenna  (SA)  method  for  operational  wind  profiling.  There  is 
virtually  no  difference  in  cost  between  the  two  types  of  systems.  However,  there  may  be 
some  significant  advantages  of  the  SA  method,  particularly  when  smaller  radars  are  being 
considered.  Since  any  Profiler  network  is  still  only  in  the  planning  stage,  now  is  the 
time  to  consider  the  various  alternatives. 

2.  Description  of  the  two  techniques 

The  Doppler  method  uses  the  same  antenna  array  for  transmitting  and 
receiving.  The  array  is  phased  in  such  a  way  that  the  beams  point  at  some  angle  off 
vertical.  The  Doppler  shift  of  the  received  signal  is  then  proportional  to  the  line-of- 
sight  velocity,  V'  in  the  right-hand-side  of  Figure  1.  Two  beams  pointing  in  different 
directions  have  to  be  used,  together  with  the  assumption  that  the  vertical  velocity  is 
negligible,  to  determine  the  two  horizontal  wind  components.  The  three-dimensional 
vector  wind  velocity  can  be  determined  uniquely  only  if  three  beam  directions  are  used. 

^ISCAT  Scientific  Association,  P.  0.  Box  705,  S-98127  Kiruna,  Sweden,  on  leave  from 
Max-Planck-Institut  fur  Aeronomie,  D— 3 1  SCatlenburg-Lir.dau,  W.  Germany 
2School  of  Electrical  Engineering,  Cornell  University,  Ithaca,  MY  14853 
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The  spaced  antenna  method  (see  Rottger,  1981a  for  a  detailed  description),  shown  in  the 
left-hand-side  of  Figure  1 ,  uses  one  transmitter  array  and  three  closely  spaced 
receiving  arrays  with  all  beams  pointing  vertically.  The  three  receiving  antennas  may 
be  separated  from  the  transmitting  antenna,  but  more  efficient  is  a  configuration  in 
which  each  third  of  the  transmitting  array  is  used  as  a  separate  receiving  array.  The 
horizontal  velocities  are  calculated  using  the  time  lags  at  which  the  cross-correlations 
of  the  signals  received  in  the  various  antennas  maximize.  The  spaced  antenna  method 
essentially  tracks  the  propagation  of  a  perturbation  in  the  refractive  index  across  the 
distance  separating  the  receiving  antennas.  The  vertical  velocity  is  calculated  from 
the  Doppler  shift. 

An  objection  to  the  spaced  antenna  method  that  is  often  heard  is  that  the 
technique  gives  incorrect  wind  measurements  when  the  medium  is  strongly  affected  by  wave 
activity.  In  that  case  the  scattering  is  primarily  from  the  refractive  index  structure 
associated  with  the  wave  fronts,  and  the  measured  velocity  is  the  phase  velocity  of  the 
waves  rather  than  the  true  wind  velocity.  However,  Briggs  (I960)  has  shown  rigorously 
that  the  Doppler  method  and  spaced  antenna  method  are,  in  fact,  equivalent.  Both 
techniques  scatter  from  the  same  variations  in  the  medium.  Therefore,  the  effect  of 

wave  structures  is  equally  a  problem  with  the  Doppler  method  and  the  SA  method.  There 
is  no  reason  to  expect  that  the  latter  will  give  poorer  results  on  that  account. 


3-OIM  VELOCITY  MEASUREMENTS  WITH  VHF-RAOAR 


Figure  1.  Schematic  of  the  spaced  antenna  and  the  Doppler  method  for  measuring  the 
three-dimensional  velocity  in  the  troposphere,  stratosphere,  and  mesosphere. 


3.  A a  experiment  to  compere  the  two  techniques 


The  SOUST-VHF-Radar  is  one  of  the  few  radars  at  present  where  a  comparison 
between  the  Doppler  technique  and  the  spaced  antenna  technique  can  be  carried  out 
directly.  The  configuration  of  the  antenna  system  is  shown  in  Figure  2.  The  large 
antenna  can  be  used  for  both  transmitting  and  receiving  in  the  Doppler  mode.  The  three 
smaller  antennas  can  also  be  used  as  spaced  receivers  in  combination  with  the  large 
transmitting  array.  In  October  1979»  an  experiment  was  carried  out  with  this  radar  to 
compare  the  two  measurement  techniques. 

The  spaced  antenna  wind  measurements  can  be  made  either  by  cross-correlating 
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Figure  2.  Antenna  systems  of  the  S0CSY-7HF- Radar  (from  Rottger  1981a).  Antenna  I  is 
used  in  the  Doppler  mode  for  transmission  and  reception.  Antenna  II  is  the  smaller 
system  of  3  x  32  four-element  Yagi  antennas  which  are  used  in  the  spaced  antenna  mode 
for  reception  when  the  large  array  of  196  lagis  is  used  for  transmission. 


the  received  power  at  the  three  receiving  antennas  or  by  cross-correlating  the  complex 
amplitudes.  The  latter  approach  retains  the  phasing  information  and  gives  better 
results  because  it  is  consistent  with  the  process  of  coherently  integrating  the  signals. 
Figure  3  shows  the  resulting  comparison.  It  is  clear  that  the  complex  amplitude 
cross-correlation  provides  data  over  a  greater  height  range  than  the  simple  power  cross¬ 
correlation.  Especially,  the  region  between  8.2-km  and  9.7-km  altitude  and  the  lower 
stratospheric  heights  should  be  compared. 

It  is  well  known  that  signals  at  7RF  have  a  strong  angular  dependence  near 
the  zenith.  Specular  reflections  or  Fresnel  scattering  associated  with  the  stable 
temperature  stratification  in  the  atmosphere  produce  enhanced  echoes  at  vertical 
incidence.  The  echo  strength  decreases  rapidly,  by  as  much  as  10  dS  (Green  and  Gage, 
1980),  at  even  a  few  degrees  off  vertical.  The  spaced  antenna  technique  makes  maximum 
use  of  this  effect  since  all  the  beams  are  pointed  vertically,  producing  a  higher 
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signal- to- noise  ratio  at  those  heights  where  the  winds  are  measured.  The  Doppler  method 
does  not  take  advantage  of  the  effect  since  the  beams  have  to  be  pointed  off  vertical  to 
make  the  wind  measurements. 

The  strong  angular  dependence  of  the  reflectivities  does  pose  a  potential 
problem  for  the  Doppler  technique  when  small  antenna  arrays  are  used.  There  are 
generally  practical  limitations  on  how  far  off  vertical  the  beams  can  be  pointed,  and, 
of  course,  the  beam  width  increases  as  the  dimension  of  the  antenna  array  decreases. 
When  the  beam  width  is  wide,  the  received  signal  is  really  a  convolution  of  the  antenna 
pattern  and  the  angular  dependence  of  the  echo  strength,  as  shown  in  Figure  4.  The 
result  is  that  the  apparent  look  angle  is  different  than  the  real  look  angle,  and  the 
horizontal  winds  tend  to  be  underestimated.  The  effect  is  critical  when  the  vertical 
pointing  direction  is  within  the  main  lobe  of  the  radiation  pattern. 


3  OCT  vm  1021- lost  GMT 


Figure  3.  Height  profiles  of  drift  speed  deduced  from  (a)  power  cross-correlation  and 
(b)  complex  amplitude  cross-correlation. 


Figure  4.  Formation  of  an  apparent  beam  direction. 
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The  effect  la  shown  in  the  experimental  results  presented  in  Figure  5. 
Figure  5a  shows  the  velocity  profile  derived  from  the  complex  crosscorrelation  of  the 
spaced  antenna  signals.  Simultaneous  aircraft  wind  measurements,  shown  by  the  open 
circles,  Indicate  good  agreement  between  the  two  independent  wind  measurements.  Panel  b 
shows  the  vertical  velocity  profile  which  could  only  be  measured  with  the  VHF  radar. 
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Figure  5.  Height  profiles  of  (a)  wind  direction  and  wind  speed  measured  with  the  spaced 
antenna  method  and  comparison  with  aircraft  winds  shown  by  the  open  circles,  (b) 
vertical  velocity,  and  ( c )  wind  speed  measured  with  the  Doppler  method  using  a  beam 
pointing  at  7°  and  (d)  3.5°  off  vertical. 

Figure  5c  gives  the  velocity  profile  derived  from  the  Doppler  method  using  a  beam 
pointed  7°  off  vertical  and  shows  very  good  agreement  with  the  spaced  antenna  method. 
However,  the  beam  width  of  the  large  array  is  5°>  and  when  the  beam  is  pointed  at  only 
3.5°  off  vertical,  there  i3  a  significant  contribution  from  the  vertical  direction.  The 
resulting  wind  profile  i3  much  poorer  since  it  has  broader  scatter  and  underestimates 
the  velocities. 

4.  Discussion 

The  data  from  the  October  1979  experiment  are  not  conclusive,  but  they  do 
indicate  that  there  may  be  certain  advantages  in  using  the  spaced  antenna  method, 
particularly  for  systems  with  small  dimensions  of  the  type  likely  to  be  used  in 
operational  wind  profiling.  Problems  associated  with  aspect  sensitivity  only  arise  when 
the  beam  width  is  large.  It  is,  in  fact,  likely  that  when  larger  systems  are  used, 
there  is  little  difference  between  the  two  techniques  since  the  beam  widths  are  small 
then.  Therefore,  the  SOWS!  comparison  may  not  fully  show  the  advantages  of  the  spaced 


antenna  method  for  small  tropospheric  sounding  systems,  and  further  comparisons  should 
be  carried  out  in  the  future. 
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Abstract.  In  February  of  1982  the  SOUSY- 
VHF-Radar,  located  near  Hannover,  West  Germany, 
was  operated  during  the  passage  of  a  warm 
front.  The  timing  of  the  radiosonde  ascent 
from  nearby  Hannover  was  such  that  the  narrow 
frontal  zone  boundary  was  immediately  over¬ 
head  at  a  height  of  approximately  7  km.  The 
data  from  this  event  has  made  it  possible  to 
make  a  detailed  comparison  of  the  radar  and 
radiosonde  reflectivity  data  that  has  not  been 
possible  before.  We  find  that  the  frontal 
boundary  thac  is  detected  by  the  radar  appears 
as  a  very  low  tropopause  height  in  the  radio¬ 
sonde  data  (approximately  3  km  less  than  the 
mean  during  this  time) . 

Introduction 

Past  studies  of  frontal  zone  dynamics  have 
relied  either  on  research  aircraft  (e.g-, 
Danielsen,  1968;  Shapiro,  1974,  1978,  1980, 

1981;  Shapiro  and  Kennedy,  1981)  or  cm-wave- 
length  Doppler  radar  observations  (e.g.. 

Browning  and  Harrold,  1970;  James  and  Browning, 
1979;  Hobbs  and  Biswas,  1979;  Hobbs  and 
Persson,  1982;  Carbone,  1982)  to  resolve  the 
small-scale  features.  Recently,  a  number  of 
frontal  passage  observations  have  been  made 
with  the  SOUSY-VHF-Radar  located  in  the  Harz 
Mountains  near  Hannover  in  West  Germany 
(Rottger,  1979,  1981a, b;  Rottger  and  Schmidt, 
1981).  The  radar  operates  at  a  wavelength  of 
5.6  m  and  measures  both  the  horizontal  and 
vertical  velocities  directly  overhead  with  a 
height  resolution  of  150  m  and  a  time  reso¬ 
lution  of  the  order  of  a  minute.  Wind  measure¬ 
ments  are  routinely  available  up  to  an  altitude 
of  almost  20  km.  The  S0USY  radar  has  a  fixed 
dipole  array  constrained  to  look  vertically  and 
can  obtain  data  either  through  the  Spaced 
Antenna  (SA)  or  Doppler  Beam  Swinging  (DBS) 
method.  The  following  data  were  obtained  with 
the  SA  mechod.  Thus,  the  radar  can  provide 
time/height  cross  sections,  but  it  cannot 
measure  horizontal  variations  directly  as  the 
aircraft  and  scanning  radar  can.  However,  the 
frontal-passage  observations  are  unique  because 
of  the  excellent  spatial  and  temporal  resolu¬ 
tion,  the  height  range  covered,  and  because  the 
radar  can  actually  locate  the  frontal  boundary 
as  a  band  of  increased  echo  strengths. 

The  reflectivity  measured  with  a  vertically 
beaming  6  m  radar  is  particularly  sensitive  to 
inversions  and  other  stratified  temperature 
structure  in  the  atmosphere  (Green  and  Gage, 
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1980;  Rastogi  and  Rottger,  1982)  .  Wavelengths 
in  the  range  from  1  to  10  cm  that  have  been 
used  for  meteorological  radars  traditionally 
are  also  sensitive  to  such  structure  to  some 
extent  (see  e.g..  Atlas  et  al.,  1966),  but  the 
effect  is  much  greater  at  frequencies  in  the 
VHF  range.  Also,  the  longer  wavelength  is  much 
less  sensitive  to  echoes  from  precipitation 
than  wavelengths  in  the  centimeter  range. 

In  February  of  1982  a  cold  front  and  a  warm 
front  passed  the  location  of  the  SOUSY  radar, 
and  data  were  taken  over  a  period  of  several 
days  both  preceding  and  following  the  frontal 
passage.  This  particular  event  was  unique 
because  the  radiosonde  ascent  from  nearby 
'lannover  occurred  so  that  the  radiosonde  passed 
directly  through  the  frontal  boundary  during 
the  ascent.  The  boundary  is  a  very  narrow  fea¬ 
ture  of  small  horizontal  extent  in  the  direc¬ 
tion  parallel  to  the  direction  of  movement  of 
the  front.  In  the  February  event,  the  trav¬ 
ersal  time  of  the  frontal  boundary  was  between 
2  and  6  hours,  dependent  on  altitude.  Typical¬ 
ly,  the  radiosonde  ascents  will  miss  the  bound¬ 
ary  since  successive  launches  are  separated  by 
12  hours.  The  timing  in  this  particular  event 
was  fortuitous,  and  the  comparisons  between  the 
radar  and  radiosonde  data  show  good  agreement. 
We  will  concentrate  on  a  comparison  of  the 
radar  reflectivity  data  and  the  reflectivity 
calculated  from  the  nearby  Hannover  radiosonde 
data. 

Synoptic  Situation 

On  February  7,  1982  a  secondary  low- 
pressure  center  formed  along  a  long  occluded 
frontal  boundary  associated  with  a  low  centered 
over  Iceland.  The  new  low  became  identifiable 
approximately  200  km  west  of  Brest,  France  at 
12Z  on  February  7  (Deutscher  Wetterdienst, 

1982a, b).  By  this  time  both  che  warm  and  cold 
front  were  identifiable  at  che  surface,  though 
the  horizontal  temperature  gradients  were  weak. 
The  300-mb  and  500-mb  upper  air  charts  showed  a 
cutoff  low  15  due  west  of  Morocco.  At  12Z  a 
short  wave  with  main  axis  running  from  north¬ 
west  to  southeast  between  40  N  and  50  N  was 
deepening  at  300-mb  as  the  trough  axis  aligned 
with  the  cutoff  low.  During  the  next  24  hours, 
che  short  wave  had  moved  approximately  20  in 
longicude  further  eastward,  and  the  wave  axis 
lay  along  the  western  coast  of  France.  By  OOZ 
on  February  9,  significant  amplification  had 
taken  place,  and  a  cutoff  high  had  formed  over 
the  Bay  of  Biscay.  The  winds  aloft  over  Cen¬ 
tral  Europe  had  changed  direction  from  essenti¬ 
ally  zonal  to  a  strong  northerly  flow. 

Figure  l  shows  the  time  height  cross 
section  of  temperature  as  measured  by  the  Han- 
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Fig.  1.  Time/height  cross-section  of  temper¬ 
ature  from  the  Hannover  radiosonde  station. 
Heavy  solid  lines  show  the  boundaries  between 
warm  and  cold  air.  The  J's  show  the  location 
of  the  jets  as  determined  from  the  radar  wind 
measurements . 

nover  radiosonde.  Regions  of  warm  and  cold  air 
have  been  labeled,  and  the  locations  of  the 
fronts  separating  the  different  air  masses  are 
indicated.  Both  fronts  are  fairly  weak,  but 
the  warm  front  is  narrower  in  horizontal  ex¬ 
tent  than  the  cold  fronc.  The  position  of  the 
two  jets,  indicated  by  "J's",  was  determined 
from  the  radar  wind  measurements  that  will  be 
presented  in  more  detail  in  a  later  paper. 

Radar  Data 

Radar  data  are  available  for  the  period  be¬ 
ginning  at  00Z  on  February  6  and  ending  at  00Z 
on  February  11.  The  radar  measures  signals 
caused  by  partial  reflections  from  turbulence 
induced  structures  in  the  refractive  index. 

One  iarge  transmitting  antenna  was  used.  To 
measure  horizontal  winds  three  separate  arrays 
were  used  to  receive  the  signals,  and  the  hori¬ 
zontal  winds  were  calculated  from  the  cross- 
correlation  between  the  echoes  received  at  each 
array.  The  vertical  winds  were  measured  dir¬ 
ectly  from  che  Doppler  shift  of  the  echoes. 

The  echo  power  is  the  third  parameter  measured 
by  che  radar,  3nd  ic  is  converted  into  relative 
radar  reflectivity  by  correcting  for  che  r^ 
dependence.  The  technique  is  described  by 
Rot  cger  ( 1981c) . 

During  che  period  discussed  here,  measure¬ 
ments  were  made  for  i  min  every  hour  on  che 
hour.  Ail  data  presented  are  an  average  value 
cor  che  corresponding  d-min  period.  This  par¬ 
ticular  sampling  scheme  was  used  to  reduce  the 
amount  of  raw  daca  thac  had  to  be  processed  co 
one  cape  per  2A  hours.  Ac  the  present  time 
only  ore inceerat ion,  buc  no  further  preproces¬ 
sing,  of  che  received  signals  is  done  on  line. 

The  height  range  in  which  there  is  a  usable 
signal  stretches  from  the  lower  limit  of  1.33 
km  MSL,  900  m  above  the  radar,  to  a  maximum  of 
approximately  20  km.  The  lower  limit  is  set  by 
che  receiver  recovery  time,  ground  clutter,  and 
atenna  near-fieid  effects.  The  upper  limit  is 
determined  by  the  decreasing  signai-to-noise 
ratio  obtained  with  the  applied  average  trans¬ 
mitter  power  of  20  KW,  the  effective  antenna 
aperture  of  2500  m“ ,  and  che  given  range  reso¬ 
lution  of  150  m. 


Reflectivity  Data 

The  reflectivities  measured  with  the  VHF 
radar  as  a  function  of  time  and  height  are 
shown  in  Figure  2.  The  increment  between  con¬ 
tours  is  3  dB,  and  the  regions  of  stronger 
echoes  have  been  shaded.  A  band  of  higher  re¬ 
flectivities  is  evident  above  11  or  12  km  in 
the  vicinity  of  the  cropopause.  Long-wave- 
length  radars  are  extremely  sensitive  to  strat¬ 
ification  in  the  temperature  structure  of  the 
atmosphere,  and  it  has  been  shown  by  Rottger 
(1980)  and  Green  and  Gage  (1980)  that  radars 
operating  in  the  VHF  regime  can  be  used  to 
determine  the  height  of  che  tropopause.  The 
radiosonde  station  located  closest  to  the  SOUSY 
site  is  Hannover,  100  km  northwest.  The  tropo¬ 
pause  heights  determined  from  the  radiosonde 
daca  at  the  standard  12-hour  intervals  have 
been  plotted  in  Figure  2  as  the  heavy  crosses. 
There  is  good  agreement  between  the  radar  and 
radiosonde  observations  taking  into  account 
chat  che  radar  sees  the  "radar"  tropopause 
which  is  a  few  hundred  meters  above  che  meteor¬ 
ological  tropopause  (Rastogi  and  Rottger, 

1982). 

At  12Z  on  February  8,  a  band  of  strong  re¬ 
flectivities  stretches  downward  from  the  tropo¬ 
pause  and  merges  with  the  region  of  strong 
echoes  below  7  km.  The  band  coincides  with  the 
location  of  the  warm  front  shown  in  Figure  1. 
The  cold  front  does  not  show  up  as  clearly,  but 
if  the  shading  threshold  is  decreased  by  6  dB, 
as  shown  by  the  cross-hatched  area  in  Figure  1, 
the  structure  associated  with  the  cold  front 
begins  to  emerge.  Presumably  the  4-min 
sampling  period  is  too  short  to  smooth  out  the 
short  term  fluctuations  of  the  weaker  signals 
from  che  cold  frontal  boundary.  The  radar  is 
able  co  detect  the  frontal  surface  for  the  same 
reason  that  it  can  detect  the  lower  edge  of  the 
tropopause  boundary.  The  temperature  stratifi¬ 
cation  associated  with  the  front  enhances  the 
echoes  at  long  wavelengths.  This  type  of 
structure  has  been  observed  consistently  with 
the  radar  during  a  number  of  ocher  warm  frontal 
passages.  In  this  particular  case  we  have  much 
more  detailed  corroboration  from  the  indepen¬ 
dent  raciosor.de  measurements,  showing  tnat  the 
structure  is  real.  Figure  2  shows  that  che 
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Fig.  2.  Time/height  cross-section  of  the 
reflectivities  measured  by  che  SOL'S': -VKF-Radar . 
The  contour  interval  is  3  dB.  Regions  of  nigh- 
est  echo  intensity  have  been  stippled.  The 
cross-hatched  areas  show  echoes  that  are  within 
b  dB  of  che  echo  strengths  in  the  stippled 
areas . 
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Fig.  3.  Logarithm  of  che  vertical  gradient  of 
the  potential  refractive  index  calculated  from 
che  Hannover  radiosonde  data. 


frontal  boundary,  i.e.,  che  boundary  between 
the  stippled  and  nonstippled  area,  is  extremely 
narrow  in  time  and  therefore,  we  deduce,  also 
in  the  horizontal  direction.  The  twelve  hour 
interval  between  successive  balloon  launches 
makes  it  difficult  to  observe  with  che  conven¬ 
tional  meteorological  observing  network.  How¬ 
ever,  in  this  case  che  boundary  passage  coin¬ 
cided  with  the  balloon  launch  at  12Z.  The 
tropopause  height  reported  by  Hannover  at  chat 
time  was  extremely  low  but  in  agreement  with 
che  radar  daca,  giving  a  strong  indication  that 
the  band  of  higher  radar  reflectivities  really 
shows  the  position  of  che  front. 

The  radar  reflectivities  are  proportional 
to  M2  with  M  given  by  (Ottersten,  1969) 


,,  ,  ,A— 0  ?  -Sine, 

-7  /  .6  x  10  —  i— ; — ; 

l  3Z 


Hr.q/jz, 
3  Ini  3z 


Here  M  is  che  mean  vertical  gradient  of  the 
potential  refractive  index,  ?  is  pressure,  T  is 
temperature,  q  is  the  specific  humidity,  and  1 
is  the  potential  temperature.  The  radiosonde 
daca  from  Hannover  were  used  to  calculate  the 
reflectivities  using  the  above  formula,  and  the 
results  are  plotted  in  Figure  3.  The  contour 
interval  is  3  dB.  The  same  features  are  evi¬ 
dent  in  both  Figures  2  and  3,  though  the 
measured  radar  reflectivities  show  more  of  the 
detailed  structure  since  che  time  resolution  is 
more  than  a  factor  of  ten  better. 


Conclusion 

The  data  presented  here  have  shown  that  the 
radar  can  detect  the  warm  frontal  boundary  even 
when  the  temperature  gradients  are  relatively 
weak.  In  general  there  is  a  greac  pocer.tial 
for  using  sensitive  radars  like  the  50USY-VHF- 
Radar  for  studies  of  synoptic  and  mesoscale 
interactions.  The  high  time  and  spatial  reso¬ 
lution  that  can  be  achieved  roucinely  with  such 
an  instrument  and  the  relatively  low  cose  of 
operation  once  the  initial  investment  for  con¬ 
struction  of  che  radar  facility  has  been 
covered  are  two  features  of  this  measurement 
technique  that  cannot  be  rivaled  by  other 


techniques.  Larsen  and  Rottger  (1982)  have 
discussed  the  applications  of  the  radars  to 
meteorological  research  in  greater  detail  . 

Acknowledgments .  The  assistance  of  staff 
members  of  the  Max-Planck-Instituc  is  grate¬ 
fully  acknowledged,  particularly  G.  Monecke . 

The  authors  are  thankful  for  the  support  of 
the  Deutscher  Wetterdiensc  in  providing  the 
surface  and  upper  air  charts,  and  che  radio¬ 
sonde  data.  MFL  was  supported  by  the  Air  Force 
Office  of  Scientific  Research  under  grant 
AFOSR-80-0020  during  the  course  of  this  study. 

References 

Atlas,  K. ,  K.  R.  Hardy,  K.  M.  Glover,  I.  Katz, 
and  T.  G.  Konrad,  Tropopause  detected  by 
radar,  Science ,  133 ,  1110-1112,  1966. 

3alsley,  3.  3.,  and  K.  S.  Gage,  On  the  use  of 
radars  for  operational  wind  profiling,  3ull . 
Amer .  Met  ■  Soc  . ,  63 ,  1009-1013,  1982. 
3rowning,  K.  A.,  and  T.  W.  Harrold,  Air  motion 
and  precipitation  growth  at  a  cold  front. 
Quart .  J.  Roy ■  Meteor ■  Soc ■ ,  96 ,  369-389, 
1970. 

Carbone,  R.  E.,  A  severe  frontal  rair.band.  Fart 
I:  Stormwide  hydrodynamic  structure,  J. 

Atmos .  Sci ■ ,  39 ,  253-279,  1982. 

Danielsen,  E.  F.,  Stratospheric-tropospheric 
exchange  based  on  radioactivity,  ozone,  and 
potential  vorticity,  2_.  Atmos  ■  Sci.,  23,  502- 
518,  1968. 

Deutscher  Vietterdienst ,  We  1 1  erkar  te  ies 

Deutschen  Wet cero i er.s :es .  Deutscher  Wetrer- 
dienst,  Seeve t reramr ,  Hamburg,  1982a. 
Deutscher  Wetterdier.st ,  Europai scher 
Vietterbericht ,  Deutscher  Wetterdiensc, 
Zentralamt,  Offenbach,  1982b. 

Green,  J.  L.,  and  X.  S.  Gage,  Observations  of 
stable  layers  in  the  troposphere  and  strato¬ 
sphere  using  VHF  radar,  Radio  Sci  . ,  15 ,  395- 
406,  1980. 

Hobbs,  P.  V.,  and  K.  R.  3iswas,  The  cellular 
structure  of  narrow  cola-frontal  rainbands, 
Quart .  _J.  Roy ■  Meteor ■  Soc  .  .  103  ,  723-727  , 
1979. 

Hobbs,  ?.  V.,  ar.d  ?.  0.  G.  Persson,  The  meso¬ 
scale  and  microscale  structure  and  organiza¬ 
tion  of  clouds  and  precipitation  in  mic- 
latitude  cyclones.  Part  V:  The  substruc¬ 
ture  of  narrow  cold-frontal  rainbands,  J. 
Atmos .  Sci.,  39 ,  230-295,  1982. 

James,  ?.  K.,  and  K.  A.  Browning,  Mesoscale 
structure  of  line  convection  at  surface  cold 
fronts.  Quart.  J.  Rov.  Meteor.  Soc.,  105, 
371-382,  197°. 

Larsen,  M.  F.,  ar.d  J.  Rottger,  VHF  and  VHF 
Doppler  radars  as  tools  for  synoptic  meteor¬ 
ology,  3ull.  Amer.  Met.  Soc.,  63,  996-1008, 
1982. 

Ottersten,  H.,  Mean  vertical  gracient  of 

potential  refractive  index  in  turbulent  mix¬ 
ing  and  radar  detection  of  OAT,  Radio  Sci., 

4,  12-7-1249,  i<369. 

Rascogi,  P.  K.,  and  J.  Rottger,  VHF  radar  ob¬ 
servations  of  coherent  reflections  in  the 
vicinity  of  the  tropopause,  J.  Atmos.  Terr. 
Phvs  ■ ,  ,  -61-469,  1982. 

Rottger,  J.,  VHF  radar  observations  of  a 

frontal  passage,  J.  Appl.  Meteor.,  13,  35-91, 
1979. 


328 


Larsen  and  Rottger:  Radar  Location  of  Frontal  Boundary 


Rottger,  J.,  Structure  and  dynamics  of  the 
stratosphere  and  mesosphere  revealed  by  VHF 
radar  investigations,  Pure  Appl.  Geophys . , 
US,  494-527,  1980. 

Rottger,  J.,  The  capabilities  of  VHF  radars  for 
meteorological  observations,  in  Nowcasting: 
Mesoscale  Observations  and  Short-Range  Pre¬ 
diction,  European  Space  Agency,  Paris,  1981a. 

Rottger,  J.,  Wind  variability  in  the  strato¬ 
sphere  deduced  from  spaced  antenna  VHF  radar 
measurements.  Preprint  Vo 1 ume ,  20th  Confer¬ 
ence  on  Radar  Meteorology,  Nov.  30-Dec.  3, 
1981,  (Boston),  AMS,  Boston,  Mass.,  1981b. 

Rottger,  J.,  Investigations  of  lower  and  middle 
atmosphere  dynamics  with  spaced  antenna 
drifts,  J.  Atmos.  Terr.  Phys . ,  43,  277-292, 
1981c. 

Rottger,  J.,  and  G.  Schmidt,  Characteristics  of 
frontal  zones  determined  from  spaced  antenna 
VHF  radar  observations.  Preprint  Volume ,  20th 
Conference  on  Radar  Meteorology ,  Nov.  30- 
Dec.  3,  1981  (Boston),  AMS,  Boston,  Mass., 
1981. 

Shapiro,  M.  A.,  A  multiple  structured  frontal 


zone-jet  stream  system  as  revealed  by  meteor¬ 
ologically  instrumented  aircraft,  Mon.  Wea. 
Rev.,  102,  244-253,  1974. 

Shapiro,  M.  A.,  Further  evidence  of  the  meso¬ 
scale  and  turbulent  structure  of  upper  level 
jet  stream-frontal  zone  systems,  Mon.  Wea. 
Rev.,  106,  1100-1111,  1978. 

Shapiro,  M.  A.,  Turbulent  mixing  within  tropo- 
pause  folds  as  a  mechanism  for  the  exchange 
of  chemical  constituents  becween  the  strato¬ 
sphere  and  troposphere,  J.  Atmos .  Sci  . ,  37 , 
994-1004,  1980. 

Shapiro,  M.  A.,  Frontogenes is  and  geostrophic- 
ally  forced  secondary  circulations  in  the 
vicinity  of  jet  stream-frontal  zone  systems, 
J_.  Atmos .  Sci . ,  38,  954-973,  1981. 

Shapiro,  M.  A.,  and  P.  J.  Kennedy,  Research 
aircraft  measurements  of  jet  stream  geo- 
strophic  and  ageostrophic  winds,  -J-  Atmos ■ 
Sci.,  38,  2642-2652,  1981. 


(Received  December  20,  1982; 
accepted  February  1,  1983.) 


Invited  review:  Proceedings  of  the  SCOSTEP/ 
URSI/MAP  Workshop  on  Technical  Aspects  of 
MST  Radar,  May  23-27,  1983,  Handbook  for 
MAP  series,  University  of  Illinois,  Urbana, 
Illinois,  1983.  . 


The  MST  Radar  Technique: 

Requirements  for  Operational  Weather  Forecasting 

Miguel  Folkmar  Larsen 
School  of  Electrical  Engineering 
Cornell  University 
Ithaca,  MY  14853 

1 .  Introduction 

Recent  interest  in  improving  our  ability  to  make  mesoscale  forecasts  has 
been  the  result  of  improving  observational  technology  and  theoretical 
understanding  of  mesoscale  motions.  There  Is  a  feeling  that  the  accuracy  of 
forecasts  for  spatial  scales  of  less  than  1000  km  and  time  scales  of  less  than 
12  hours  can  be  improved  significantly  if  our  resources  are  applied  to  the 
problem  in  an  intensive  effort  over  the  next  decade.  Since  the  most  dangerous 
and  damaging  types  of  weather  occur  at  these  scales,  there  are  major  advantages 
to  be  gained  if  such  a  program  is  successful.  The  interest  in  improving  short 
term  forecasting  is  evident  in  the  series  of  papers  resulting  from  the 
International  IAMAP  Symposium  on  Howcasting  [3rowning,  1982],  the  paper  on  the 
long-term  goals  of  NCAA  as  presented  by  Schmidt  [1583]  at  the  Fifth  Symposium 
on  Meteorological  Observations  and  Instrumentation,  and  the  two  handbooks 
published  by  the  UCAR  committee  on  the  National  STORM  (Stormscale  Operational 
and  Research  Meteorology)  Program  [UCAR,  1982;  Anthes,  1983].  The  conclusion 
of  all  of  these  publications  is  that  the  technology  at  the  present  time  is 
sufficiently  developed,  both  in  terms  of  new  observing  systems  and  the 
computing  power  to  handle  the  observations,  to  warrant  an  intensive  effort  to 
improve  stormscale  forecasting. 

The  National  STORM  Program  handbooks  [UCAR,  1982,  Anthes,  1983]  are 
excellent  source  books  on  the  general  problem  of  short  term  forecasting.  The 
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first  handbook  summarizes  the  need  for  this  type  of  forecasting  and  the 
framework  for  achieving  the  desired  improvements.  The  second  handbook  gives  a 
detailed  assessment  of  our  capabilities  and  understanding  at  the  present  and 
the  areas  of  research  most  in  need  of  attention.  Questions  that  will  have  to 
be  addressed  deal  with  the  type  of  observing  system  or  systems  that  will  be 
used  to  generate  data  compatible  with  the  scales  of  motion  being  discussed  and 
the  way  the  data  will  be  used  operationally.  Hooke  [1983]  summarized  the 
situation  as  follows: 

Within  the  next  several  years,  operational  meteorology  will  be 
facing  some  crucial  decisions  regarding  measurement  systems.  The 
most  obvious  and  important  example  is  the  determination  of  winds.  A 
variety  of  remote  sensing  techniques  have  been  used  experimentally  to 
measure  wind  direction  and  speed  in  clear  air.  There  is  substantial 
need  for  an  integrated  effort  to  determine  which  of  these  techniques 
is  most  promising  by  testing  alternative  methods  in  the  field.  This 
will  be  a  necessary  prelude  to  the  procurement  and  deployment  of  a 
next-generation  wind  sensing  system. 

In  thi3  review  I  will  attempt  to  provide  an  assessment  of  the  extent  to  which 
the  so-called  MST  radar  technique  fulfills  the  requirements  for  an  operational 
mesoscale  observing  network,  and  I  will  delineate  the  extent  to  which 
improvements  in  various  types  of  forecasting  could  be  expected  if  such  a 
network  is  put  into  operation. 

2.  Characteristics  and  capabilities  of  the  radar  technique 
The  MST  radar  measures  the  vertical  profile  of  horizontal  winds  based 


either  on  the  Doppler  shift  of  the  signal  backscattered  from  turbulent 
variations  in  the  refractive  index  or  from  the  cross-correlation  between  the 
signals  received  at  a  set  of  three  spaced  antennas.  The  technique  has  been 
described  in  detail  by  Woodman  and  Guillen  [1974],  Gage  and  Balsley  [1973], 
Balsley  and  Gage  [1990],  and  Rottger  [1980]  among  others.  Details  and 
advantages  of  various  specific  measurement  techniques  are  discussed  in  much 
greater  detail  in  other  papers  in  this  volume. 

I  will  concentrate  on  VHF  radars  that  use  a  fixed  dipole  array  and  are 
capable  of  measuring  the  vertical  profile  of  the  horizontal  winds  as  a  minimum. 
Some  MST  radars  operating  at  shorter  wavelengths  have  fully  steerable  dishes. 
An  example  is  the  23  cm  radar  formerly  located  at  Chatanika,  Alaska,  and  now 
located  in  Stfndre  Strsdmfjord,  Greenland.  I  will  not  consider  such  radars  here 
since  they  are  in  the  same  category  as  the  existing  weather  radar  network  if 
upgraded  to  provide  Doppler  capability. 

The  radar's  ability  to  provide  wind  profiles  is  something  that  is 
duplicated  by  the  standard  rawinsonde.  However,  the  radar  wind  measurements 
have  many  advantages,  particularly  if  the  mesoscale  is  of  interest.  The  radar 
wind  profiles  can  be  measured  as  often  as  desired  down  to  the  imposed  equipment 
limitation  of  a  few  minutes.  There  are  no  expendables  involved  except  for  the 
power  used  to  run  the  radar  equipment.  Also,  the  radar  measures  the  wind 
profile  immediately  above  the  radar.  During  periods  of  high  winds,  the 
rawinsonde  can  drift  as  much  as  100  km  during  its  ascent.  Such  an  error  can  be 
a  ma.jor  one  if  the  scales  of  motion  of  interest  are  of  the  order  of  a  few 
hundred  kilometers.  Finally,  a  large  part  of  the  cost  of  running  a  radiosonde 
station  has  always  been  the  salaries  of  the  personnel.  Since  the  radar 
measurements  are  easily  automated,  adding  VHF  Doppler  radars  to  the  network 
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would  onl7  increase  the  required  ertra  manpower  by  a  small  fraction  of  the 
cumber  of  personnel  needed  if  an  equal  number  of  new  radiosonde  stations  were 
established.  The  success  of  such  an  approach  is  already  evidenced  by  the  Poker 
Flat  MST  Hadar  which  has  been  operating  unattended  since  1979. 

However,  the  comparison  between  the  radiosonde  and  the  racar  is  not 
really  the  most  appropriate.  The  possibility  of  establishing  a  meso3cale 
radiosonde  observing  network  has  sever  been  considered  very  seriously  since  the 
cost  is  prohibitive.  Most  likely  the  competition  will  come  from  satellites  or 
some  other  ground-based  remote  sensing  system  such  as  sodars  or  lidars.  The 
primary  advantage  of  the  radar  over  the  acoustic  sounder  is  the  height  range 
covered  and  the  insensitivity  of  the  former  to  various  forms  of  noise.  A  very 
complete  review  of  the  capabilities  of  the  acoustic  sounding  technique  is  given 
by  Brown  and  Hall  [1978].  The  lidar  is  competitive  with  the  radar  in  terms  of 
the  height  range  covered  and  the  time  and,  certainly,  the  height  resolution, 
but  the  lidar  is  seriously  hampered  when  there  is  precipitation  or  fog  or 
simply  when  it  is  overcast  [Strauch  and  Cohen,  1972].  The  lidar  provides 
information  on  atmospheric  density  and  humidity,  but  the  radar  provides  data  on 
the  height  of  the  tropopause  and  other  temperature  structure  such  a3  inversions 
and  fronts.  Further  comparison  between  the  advantages  of  the  two  techniques  is 
needed,  although  the  main  drawback  of  the  lidar  for  operational  applications 
appears  to  be  the  limitations  of  the  technique  when  used  in  the  presence  of 
clouds. 

The  satellite  has  a  number  of  advantages  over  the  racar.  Primarily  it 
provides  vast  areal  coverage  in  relatively  short  periods  of  tine.  Since  a  very 
large  source  of  error  for  large  scale  forecasting  is  due  to  a  lack  of  data  in 
sparsely  populated  regions,  underdeveloped  areas,  and  the  oceans,  satellite 
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measurements  appear  to  be  the  beat  hope  for  improving  large-scale  observations. 
The  radars  do  not  immediately  offer  any  hope  for  providing  data  in  data  sparse 
regions.  The  height  and  time  resolution,  and  the  height  range  covered  by  the 
radar  is  significantly  better  than  that  of  the  satellite,  but  the  lack  of 
spatial  coverage  of  the  radar  measurement  is  still  one  of  the  limitations  of 
the  technique. 

So  far,  I  have  compared  the  advantages  and  disadvantages  of  the  various 
techniques.  It  is  clear  that  the  radar  by  itself  cannot  replace  the  other 
types  of  measurements,  but  it  can  provide  a  relatively  inexpensive  solution  to 
the  problem  of  upgrading  the  observing  network  to  provide  aesoscale  resolution. 
Meteorological  radars  already  exist  and  provide  important  information  for  the 
forecaster.  Half-hourly  satellite  photos  are  an  important  input  to  the 
forecast  process  and  satellite  measurements  of  temperature  and  winds  are 
providing  data  for  areas  where  no  information  could  be  obtained  otherwise, 
albeit  with  less  than  desirable  height  resolution.  However,  the  VHF  radars  can 
fill  the  gaps  that  exist  in  the  present  radiosonde  network.  It  has  been  shown 
by  Hoke  and  Anthes  [19763,  Daley  and  Puri  [1980],  and  Daley  [1980]  among  others 
that  when  small  scales  of  motion  are  of  interest,  the  wind  information  is  more 
important  than  pressure  and  temperature  information.  Therefore,  the  fact  that 
the  radar  only  measures  the  winds,  and  not  also  pressure  and  temperature, 
should  not  necessarily  be  viewed  as  a  compromise  of  the  technique  when  it  is 
used  for  mescscale  observations.  I  will  discuss  this  point  in  more  detail  in 
Section  5 . 

3.  Requirements  for  an  operational  system 
Although  the  topic  of  the  workshop  is  the  MST 
(Mesosphere-Stratosphere- Troposphere)  radar,  the  ability  to  measure  winds  in 
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the  troposphere,  stratosphere,  and  mesosphere  is  not  required  of  the  systems 
used  for  xesoscale  weather  forecasting.  A  simple  "Model  T"  Radar  would  be 
sufficient.  With  thanks  to  Henry  Ford,  such  a  radar  should  be  simple  to  mass 
produce,  low  cost,  and  dependable  in  operation.  The  reductions  in  antenna  size 
and  transmitter  power  gained  by  reducing  the  design  specifications  of  the 
system  could  reduce  the  cost  of  the  system  by  as  much  as  a  factor  of  ten  when 
compared  to  an  MST  system.  That  is  a  crude  estimate  but  probably  not 
unrealistic. 

Discussions  about  the  applicability  of  the  radar  wind  measurements  to 
operational  weather  forecasting  have  generally  centered  on  applications  to 
large-scale  forecasting  with  improved  spatial  resolution  so  that  mesoscale 
phenomena  can  be  resolved.  Of  course,  it  is  important  that  this  aspect  of  the 
radar  technique  should  be  discussed,  and  that  may  be  how  the  systems  are 
ultimately  applied.  However,  there  are  intermediate  applications  for  such 
systems  that  would  allow  the  radars  to  be  phased  into  the  large-scale  observing 
network  with  possible  significant  forecast  improvements  at  each  intermediate 
step.  There  are  a  number  of  very  localized  phenomena  that  lead  to  severe 
weather  that  a  small  network  of  VKF  radars  could  be  useful  in  forecasting. 

Applications  that  come  to  mind  include  the  following.  A  network  of 
between  3  and  5  radars  distributed  around  the  Great  Lakes  could  be  used  for 
operational  forecasting  of  lake-effect  snows.  Typically,  the  snows  are 
generated  by  directional  changes  in  the  mesoscale  flew  patterns  and  eventually 
mesoscale  circulations  develop  in  response  to  the  heating  effects  of  the  lakes 
[Zipser,  19833.  It  should  be  possible  to  detect  these  small  scale  circulation 
changes  with  a  network  of  radars.  More  study  would  be  needed  to  determine  the 
real  usefulness  of  such  a  system  for  forecasting  this  very  specific  local 
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phenomenon. 

A  second  application  would  involve  use  of  a  system  of  radars  as  a 
forecast  tool  for  the  severe  Colorado  wind  storms  that  occur  everyyyear  in  the 
lee  of  the  Rocky  Mountains  [Lilly  and  Zipser,  1972;  Klemp  and  Lilly,  1975]. 
Since  a  cloudy  or  precipitating  atmosphere  is  not  necessarily  associated  with 
this  phenomena,  a  scanning  weather  radar,  even  one  with  Doppler  capability,  is 
not  particularly  useful  for  forecasting  this  type  of  event.  Actually,  the 
necessary  network  may  already  exist  in  the  form  of  the  PROFS  (Prototype 
Regional  Observing  and  Forecasting  System)  network  of  wind  profilers  operated 
by  the  Wave  Propagation  Laboratory  of  NOAA  [Strauch  et  al. ,  1982], 

Another  application  of  a  small-scale  system  of  Model  T  radars  is  for 
studies  of  the  sea  breeze  in  the  Florida  peninsula.  The  sea  breeze  develops  in 
response  to  the  diurnal  heating  cycle  and  the  temperature  differences  between 
land  and  sea.  The  vertical  circulation  that  develops  acts  as  a  modulator  of 
the  convective  activity  over  the  land  and  over  the  ocean  [Lhermitte  and  Gilet, 
1975;  Atkinson,  1981].  The  resulting  thunderstorms  can  be  very  severe  and  may 
involve  large  shears,  heavy  rainfall  rates,  hail,  and  turbulence  that  can  be  a 
hazard  to  local  aviation.  A  scanning  weather  radar  can  be  used  to  detect  the 
developing  cells,  but  since  the  lifetime  of  a  single  cell  is  from  30  to  45 
minuses,  only  a  3hort  term  warning  can  be  issued.  Clear  air  wind  measurements 
may  be  capable  of  detecting  the  buildup  of  the  conditions  leading  to  intense 
convection. 

The  possible  applications  of  the  radar  systems  for  local  forecasting  just 
named  are  only  a  few  of  the  possibilities.  The  important  point  is  that  a  3nail 
cumber  of  the  radar  systems  can  be  installed  to  provide  improved  forecasting  of 
specific  local  phenomena.  Confidence  in  the  systems  and  operational  experience 


3 


would  So  gained  before  asking  a  comni  taent  to  use  the  radars  on  a  networkwide 
basis. 

Another  area  in  which  the  radars  can  be  applied  for  forecasting  purposes 
relates  to  pollutant  dispersion.  Most  of  the  aodel3  used  to  estiaate 
dispersion  of  pollutants  are  based  on  the  Gaussian  plume  aodels  [Hanna  et  al. , 
1982].  The  two  major  input  parameters  are  am  atmospheric  stability  index 
derived  from  the  radiosonde  temperature  profile  and  the  wind  at  the  height  of 
the  center  of  the  plume.  A  major  problem  is  the  significant  diurnal  variation 
in  the  wind3  that  cannot  be  resolved  by  the  rawinsonde  measurements  made  once 
every  12  hours  [Draxler,  19831.  Only  a  small  system  that  could  measure  the 
winds  up  to  a  height  of  a  few  kilometers  would  be  needed  to  improve  the  wind 
information  data  base  significantly.  Other  locations  that  could  benefit  from 
such  a  system  would  be  airports  where  the  primary  hazard  is  from  clear  air 
turbulence  and  downdrafts  that  affect  aircraft  during  takeoff  and  landings. 
Again  only  measurements  within  the  boundary  layer  or  a  little  higher  would  be 
required.  There  are  no  doubt  other  possibilities. 

4.  Previous  work 

Very  little  work  has  been  done  to  date  dealing  directly  with 
applicability  of  the  radars  to  the  forecasting  problem.  Balsley  and  Gage 
C 1 962  3  have  discussed  considerations  for  implementation  of  an  operational  radar 
system  including  antenna  size  needed,  most  favorable  frequency  ranges,  and  the 
type  of  power  needed.  Carl3cn  and  Sundararaoan  [1932]  have  made  a  preliminary 
case  showing  that  a  few  percent  of  the  annual  fuel  consumption  of  the  airlines 
could  be  saved  if  data  from  a  wind  measuring  radar  network  was  available  for 
flight  planning  purposes.  They  indicate  that  a  detailed  3tudy  has  to  be  made 
to  determine  if  their  rough  calculations  of  potential  savings  are  correct. 


However,  the  eatings  that  they  envision  would  be  enough  to  justify  the  cost  of 
a  radar  network  within  the  first  year. 

Fukao  et  al.  [1982]  have  made  a  detailed  comparison  of  rawicsonde  data 
from  San  Juan,  Puerto  Hico  and  wind  profiles  measured  with  the  Arecibo  430  MHz 
radar.  Twenty-aim  separate  days  from  August  and  September  of  1977  were 
involved  in  the  study.  The  comparison  indicated  a  difference  of  4.9  m/s  in  the 
upper  troposphere  and  a  difference  of  3*3  n/s  in  the  lower  stratosphere.  The 
difference  in  the  lower  stratosphere  could  be  explained  by  the  experimental 
error  in  the  rawicsonde  measurement,  but  the  larger  difference  in  the  upper 
troposphere  was  apparently  due  to  spatial  variations  over  the  30  km  distance 
separating  the  two  sites. 

Larsen  [19831  investigated  the  effect  of  high  frequency  meteorological 
noise  on  the  representativeness  of  the  radar  wind  data.  If  one  Is  interested 
in  using  the  data  for  input  to  a  numerical  model  that  can  resolve  synoptic  and 
mesoscale  motions  down  to  scales  of  a  few  hundred  kilometers,  any  motions  with 
smaller  scales  are  effectively  just  noise  or  an  error  in  the  measurement.  The 
radar  wind  measurements  from  the  Poker  Flat  MST  radar  were  compared  to  the 
rawicsonde  measurements  from  Fairbanks,  Alaska,  and  both  were  compared  to  the 
geostrophic  wind  calculated  by  applying  an  objective  analysis  scheme  to  the 
standard  radiosonde  data  from  five  nearby  stations.  3y  averaging  the  high  time 
resolution  radar  data  over  intervals  of  several  hours,  the  high  frequency 
oscillations  could  be  filtered  out.  The  comparison  showed  that  the  radar  and 
rawicsonde  data  agreed  to  within  2-3  e/s  when  the  radar  data  was  averaged  over 
12  hours  or  more,  and  that  the  difference  between  the  radar  winds  and  the 
geostrophic  winds  was  about  the  sane  as  the  difference  between  the  rawicsonde 
winds  and  the  geostrophic  wind.  The  two  independent  wind  measurements  were 
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most  similar ,  and  both  differed  from  the  geoatrophic  wind  by  1. 5-2.0  m/s  more 
than  they  differed  from  each  other.  The  results  indicate  that  it  is  crucial 
that  the  data  should  be  averaged  in  some  way  if  it  is  used  as  input  to  a 
numerical  model.  The  errors  decrease  as  the  averaging  interval  is  increased. 
Therefore,  the  acceptable  error  for  a  given  model  will  have  to  be  determined  in 
order  to  know  how  to  process  the  radar  data. 

Larsen  and  Hottger  [1983]  have  shown  that  a  VEF  radar  can  detect  the 
location  of  frontal  boundary  surfaces  as  enhancements  of  the  radar 
reflectivity.  When  coupled  to  the  measurements  of  the  horizontal  and  vertical 
wind  components,  such  information  would  be  of  value  in  forecasting  the 
development  and  position  of  fronts  since  such  small-scale  features  are  not 
resolved  by  the  synoptic  observing  system  or  in  operational  numerical  forecast 
models. 

Green  et  al.  [1978]  showed  the  ability  of  the  radar  to  detect  changes  in 
the  jet  stream  height  and  intensity  in  real  time,  as  well  as  gravity  wave 
activity  and  vertical  velocities  associated  with  the  jet  stream.  Turbulence 
intensity  and  location  can  also  be  determined  along  with  the  rate  of  turbulent 
dissipation  [Gage  et  al.,  1978].  The  results  are  intriguing,  and  more  work 
needs  to  be  done  in  this  area.  One  very  interesting  aspect  of  their  study  is 
that  radar  measurements  may  be  used  in  the  future  to  provide  real  time  inputs 
for  the  parameterization  schemes  used  in  numerical  models.  All  processes  in  a 
numerical  model  with  spatial  scales  smaller  than  the  models  grid  spacing  are 
parameterized.  Such  things  a3  the  vertical  fluxes  of  heat  and  momentum  due  to 
convection  and  the  loss  of  energy  to  subgrid  scales  are  included  in  the 
paraneterizations.  The  schemes  used  usually  depend  in  seme  way  on  the  physical 
quantities  calculated  by  the  model,  but  the  high  tine  resolution  measurements 
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of  the  radar,  along  with  the  vertical  velocity  measurements  [e.g. ,  Ecklund  et 
al. ,  1982],  ma7  provide  valuable  information  on  the  actual  magnitude  of  these 
various  quantities  as  a  function  of  time  and  location. 

The  most  detailed  investigation  of  the  applicability  of  the  radar 
technique  to  mesoscale  forecasting  is  being  undertaken  by  the  Wave  Propagation 
Laboratory  of  NCAA  using  their  Profiler  system  [Strauch,  1981;  Straueh  et  al. , 
1982].  The  Profiler  uses  several  microwave  radiometers  to  measure  the 
temperature  and  humidity  profile.  A  set  of  three  Doppler  radar3  distributed  in 
a  triangular  network  around  Boulder,  Colorado  is  used  to  resolve  mesoscale 
features  in  the  local  winds.  The  system  i3  dedicated  to  the  problem  of 
mesoscale  forecasting.  Except  for  results  dealing  with  the  measurement 
capabilities  of  the  system  and  the  corresponding  accuracies,  few  results  are 
available  to  date.  However,  the  Profiler  should  provide  a  good  assessment  of 
potential  improvements  in  mesoscale  forecasting  that  can  be  realised  with  such 
a  system. 

5.  Optimal  specification  of  initial  data 

Much  of  the  discussion  about  using  the  MST  systems  operationally  has 
focused  on  the  WPL  Profiler  which  uses  a  clear  air  radar  for  wind  profiling. 
However,  the  system  was  designed  to  test  the  possibility  of  replacing  the 
standard  National  Weather  Service  radiosonde.  To  this  end,  microwave 
radicmeter3  are  used  to  provide  profiles  of  the  temperature  and  humidity, 
though  to  date  the  achievable  height  resolution  has  been  less  than  that  of  the 
radiosonde.  It  remains  to  be  seen  whether  the  difference  is  significant. 

There  is  little  doubt  that  the  thermodynamic  and  moisture  information  is 
valuable,  but  in  this  case  we  should  question  whether  the  extra  cost  would  be 
Justified  on  a  networkwide  basis.  The  radiometers  are  by  far  the  most 
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expensive  part  of  the  system.  The  cost  ratio  between  the  Model  T  radar 
discussed  here  and  the  microwave  system  may  be  as  great  as  1  to  5  for  an 
operational  system.  The  estimated  cost  of  a  Profiler  system  has  been  given  as 
approximately  $500,000  [Decker,  19831.  Carlson  and  Sundararaman  [1982] 
estimate  that  a  tropospheric  wind  profiling  system,  without  the  radiometers, 
could  be  built  for  105  of  that  sum  or  less. 

Will  the  wind  information  be  useful  in  and  of  itself?  The  answer 
apparently  is  yes.  Studies  by  Rutherford  and  Asselin  [19721,  Williamson  and 
Dickinson  [1972],  Hoke  and  Anthes  [1976],  Daley  and  Puri  [1980],  Daley  [1980], 
and  Bube  and  Ghil  [1982]  indicate  that  at  small  scales  the  wind  information  is 
by  far  the  most  important.  Small  scales  in  thi3  case  are  defined  by  a 
horizontal  scale  related  to  the  Rossby  radius  of  deformation.  The  various 
studies  have  used  either  analytic  models  or  numerical  models  that  characterize 
the  dynamics  of  the  atmosphere  and  considered  the  problem  of  how  new 
information  is  assimilated  into  the  model.  The  winds  and  the  thermodynamic 
variables  have  been  given  as  initial  conditions  singly  and  in  combination. 
Updating  the  calculated  fields  with  the  pressure  fields  would  eventually  cause 
an  adjustment  to  the  input  data,  but  the  adjustment  time  would  be  long. 
However,  the  wind  information  is  readily  absorbed  with  a  minimum  of  wave  noise 
being  generated.  The  pressure  field  then  adjusts  to  reflect  the  changes  in  the 
wind  field  so  that  a  complete  set  of  information  relevant  to  the  small-scale 
motions  is  obtained  even  though  only  the  wind  fields  are  used  as  input.  There 
seems  to  be  agreement  that  wind  information  is  mere  important  for  scales  less 
than  10C0-2CC0  km.  At  large  scales  the  reverse  is  true,  and  the  thermodynamic 
information  is  more  important. 

The  physical  reason  for  thi3  effect  is  associated  with  the  process  known 
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aa  gecstrophic  adjustment  and  first  described  by  Hossby  [1938].  The  crucial 
parameter  in  the  problem  is  the  Hossby  radius  of  deformation  given  by  the  ratio 
of  the  speed  of  sound,  or  the  propagation  velocity  of  gravity  waves  of  the 
appropriate  scale,  to  the  Coriolis  parameter  [Bluaen,  1972],  At  scales  much 
larger  than  the  Hossby  radius,  the  wind  field  will  adjust  to  balance  a 
perturbed  pressure  field.  For  scales  smaller  than  the  Hossby  radius,  the 
pressure  field  adjusts  to  balance  the  wind  field.  Therefore,  the  smaller  the 
scale  of  motion  that  is  of  interest,  the  more  useful  is  the  wind  information. 
Since  the  mesoscale  lies  in  the  range  of  spatial  scales  where  the  wind 
information  is  most  useful,  there  may  be  a  significant  improvement  in  our 
ability  to  forecast  for  this  scale  even  if  only  the  wind  fields  are  measured. 
Shapiro  et  al.  [1983]  also  discussed  the  importance  of  wind  information  at 
small  scales,  and  they  indicate  that  improved  spatial  resolution  may  be 
obtained  from  the  measurements  of  a  single  radar  if  the  data  are  used  to 
extrapolate  quantities  along  air  parcel  trajectories. 

6.  Conclusion 

This  article  was  meant  to  suggest  a  number  of  possible  applications  for 
the  MST  radar  technique  in  operational  weather  forecasting.  The  results  to 
date  have  shown  that  there  is  great  promise  for  the  technique  as  part  of  the 
standard  observing  network  once  it  is  resolved  to  increase  the  resolution  to 
include  mesoscale  motions.  In  reality,  very  little  work  has  been  done  to  date 
relating  to  the  forecast  improvements  that  actually  can  be  achieved  with  an 
operational  wind  profiling  system.  Forecast  improvements  due  to  a  neowork  of 
radars  should  be  investigaoed  both  theoretically  and  experimentally .  The 
possibility  of  using  the  radar  data  to  refine  seme  of  the  parameterization 
schemes  used  in  forecast  modeling  should  also  be  investigated.  Finally,  the 


I 
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applications  of  the  radar  systems  to  forecasting  area-specific  phenomena  should 

be  examined  in  the  near  future. 
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ABSTRACT 

Partial  and  total  aliasing  problems  in  spectral  analysis  of  the  vertical  velocity  measurement  by  single 
VHF  radar  beam  are  analyzed,  and  an  iteration  technique  for  the  aliasing  correction  is  presented.  Also,  a 
systematic  criterion  for  the  total  aliasing  correction  of  vertical  velocity  measurement  by  three  spatially 
separated  VHF  radar  beams  is  presented  and  discussed. 

Key  words:  VHF  radar  technique,  aliasing  probrem  and  its  correction 


Introduction 

Spectral  analysis  and  autocorrelation  analysis  are 
two  major  techniques  which  are  used  in  measuring  the 
vertical  velocity  by  radar.  The  velocity  u  of  the  scattering 
target  in  atmosphere  along  the  radar  beam  direction  is 
linearly  related  to  the  Doppler  frequency  /  by  v  *  y  /, 
where  \  is  the  radar  wavelength.  The  mean  velocity  uof 
the  atmosphere  layer  within  a  certain  time  period  is  natu¬ 
rally  related  to  the  mean  frequency  of  the  Doppler  spec¬ 
trum  of  the  radar  echo  signals  by  u  ■  y  /.  The  mean  fre¬ 
quencies  were  obtained  either  directly  from  the  Doppler 
spectrum  which  was  the  finite  Fourier  transform  of  the 
radar  echo  signals,  or  from  the  phase  angle  of  the  complex 
autocorrelation  function  at  the  first  time-lag,  while  the 
aliasing  problems  have  not  been  paid  enough  attention  to. 
Careful  study  of  the  periodic  nature  of  Doppler  spectrum 
led  us  to  conclude  that  the  aliasing  problems  are  serious 
enough  to  cause  grave  data  error  which  is  harmful  to  the 


phenomena  under  investigation.  In  principle,  the  aliasing 
can  be  corrected  using  a  procedure  developed  by  Sato  and 
Woodman  [  1  ] ,  however,  an  alternative  procedure  is  pre¬ 
sented  in  this  paper. 

Aliasing  in  Spectral  Analysis  and  its 
Correction 

There  exist  two  types  of  aliasing  in  the  Doppler 
spectrum,  one  is  the  nonuniform  amplification  introduced 
by  the  preintegration  of  radar  echo  signals  during  data 
acquisition  process,  the  other  one  is  the  aliasing  of  the 
frequencies  of  the  spectrum  resulted  from  the  periodic 
nature  of  finite  Fourier  transform. 

The  nonuniform  amplification  problem  can  be 
understood  as  follows.  The  output  complex  data  Z{tk ) 
at  time  :k  after  on-line  preintegration  of  n  consecutive 
signals  of  a  plane  wave  of  frequency  /,  e‘w/t,  is  given  by, 
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-g(/)exp  [i2ir/rfc.i].  (I) 


where 


?(/) 


_1_  1-exp  (ilirUt) 
n  i-exp(i2Tr/ar/«) 


for  /  =  0 
for  /^O. 


(2) 


So  the  power  of  the  original  plane  wave  with  frequency  / 
is  amplified  by  a  preintegration-amplification  factor  H(f) 
given  by 


H(f)  *  lg(/)lJ 


1  l-cos(2irfAt) 
n1  l-cos(2nfAt/n) 


(3) 


1.0 


So.5 


-0.3 


0.3 


0 

f(HZ) 

Fig.  I.  Preintegrate  amplification  factor  H(f)  at  different  fre¬ 
quencies. 


The  number  of  radar  echo  signals  n  and  time  interval  At 
of  preintegration  are  preset  by  the  experiment.  In  our 
example  data  which  will  be  described  later  in  this  section, 
n* 256  and  Af*  1.03  seconds,  the  dependence  of  H{f) 
on  the  Doppler  frequency  /  as  shown  in  Fig.  1,  clearly 
shows  that  the  non-uniformity  of  the  preintegration- 
amplification  factor  is  significant. 

Conventionally,  the  so-called  non-uniform  amplifi¬ 
cation  problem  is  really  associated  with  the  fact  that  any 
filtering  done  on  the  data  will  affect  the  data  and  to  re¬ 
store  the  original  property  of  the  data  the  transfer  func¬ 
tion  of  the  filter  should  be  used.  Coherent  integration 
amounts  to  a  low-pass  filter,  and  Fig.  1  is  indeed  the  mag¬ 
nitude  of  the  transfer  function  of  that  filter.  In  practice, 
the  coherent  integration  time  had  better  be  chosen  such 
that  it  will  not  gravely  affect  the  DoppleT  spectrum. 

The  aliasing  of  some  frequencies  of  the  Doppler 
spectrum  can  be  explained  by  the  following  analysis  of 
the  periodic  nature  of  finite  Fourier  transformation. 


Z(tk )  =  -Lr%  Z(fn)e nmklN 
yJN  ',:a0 

Z{fn)  =  ~/l0Z«k)e'7'mklN 


r*  *  KAt 

fn  ■ 


n 

NAt 


(4) 


k,  /I.  «  0,  1,2 . N- 1. 

It  can  be  easily  seen  from  Eq.  (4)  that, 

£(/„  ♦  NAf)  »  Z(/„) 
and  so  is  the  Doppier  spectrum, 

D{fn  +  NAD  *  D(fn)  (5) 

because  the  Doppler  spectrum  is  given  by 


*>(/„)»  1 

The  power  spectrum  of  the  radar  echo  scattered  from  the 
clear  atmosphere  turbulence  tends  to  be  Gaussian  which 
is  symmetric  with  respect  to  its  mean  /  as  expressed  in 
Fig.  2a.  However,  the  finite  Fourier  transform  always 
presents  the  spectrum  in  the  frequency  domain  between 
-  j  A f  and  (  j  -  1)  Af  with  its  center  at  zero  frequency, 
and  that  results  in  the  aliasing  of  some  frequencies  by 
either  NAfor  -A'A/due  to  the  periodic  nature  of  Doppler 
spectrum  as  expressed  by  Eq.  (5).  This  kind  of  aliased 
spectrum  is  demonstrated  in  Fig.  2b.  For  the  convenience 
of  further  discussion,  let’s  call  this  type  of  aliasing  as  par¬ 
tial  aliasing  of  the  Doppler  spectrum. 


Fig.  1  (a)  Original  power  spectrum. 

fb)  Partially  aliased  power  spectrum  of  fa). 
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Practically,  the  partial  aliasing  of  the  Doppler  spec¬ 
trum  as  shown  in  Fig.  2b  can  be  reconstructed  to  its  real 
form  as  shown  in  Fig.  2a  by  an  iteration  technique.  The 
iteration  process  is  as  follows,  first,  calculate  the  mean 
frequency  f  of  the  power  spectrum  obtained  directly 
from  finite  Fourier  transform,  then  shift  the  frequency 
domain  of  the  spectrum  from  [-  y  A /  (y  -1)A/]  to 
[/-  y  A/  /"+  (y  - 1)  A/]  by  applying  Eq.  (5),  followed 
by  recalculating  f  and  shifting  the  spectrum  domain  by 
putting  the  center  of  the  spectrum  at  the’newly  obtained 
mean  frequency.  Repeat  this  process  again  and  again  until 
the  mean  frequency  /"converges  to  /*. 

The  following  examples  are  from  the  result  of  our 
analysis  of  the  sample  of  SOUSY  data,  which  were  taken 
with  SOUSY-VHF  radar  by  Max-PIanck-Institute  fur 
Aeronomie.  The  radar  is  located  in  the  Harz  Mountain 
in  Northern  Germany.  Pulse  modulation,  choerent  de¬ 
tection  and  on-line  pre-integration  were  applied,  and  the 
data  were  stored  on  magnetic  tape  for  off-line  evaluation. 
The  radar  operated  at  a  frequency  of  53  MHz  and  the 
beam  was  pointed  upward  vertically.  256  consecutive 
radar  echo  signals  were  on-line  preintegrated  with  a  time 
interval  of  1.03  seconds  to  give  one  record,  each  record 
includes  105  complex  data  for  105  different  height  ranges 
between  0.9  kra  and  16.65  km,  with  height  resolution 
0.15  km.  There  are  21444  records  in  our  sample  tape. 

To  obtain  the  mean  vertical  velocity  of  each  height 
range  within  certain  time  period  by  spectral  analysis  tech¬ 
nique,  32  consecutive  complex  data  with  a  time  period 
of  32.96  seconds  (=32x1.03)  for  each  height  range  were 
analyzed  by  FFT  to  obtain  a  power  spectrum.  Then,  8 
consecutive  power  spectra  were  averaged  to  give  a  smooth¬ 
er  Doppler  spectrum,  82  such  Doppler  spectra  were  ob¬ 
tained  for  each  height  range,  and  each  Doppler  spectrum 
covers  a  time  interval  of  4.39  minutes  (=32.96  x  8  seconds). 
There  were  a  total  number  of  8610  (=82x105)  Doppler 
spectra,  among  them  4  sample  spectra  are  demonstrated 
in  Figs.  3  (a,b,c,d),  where  partial  aliasing  can  be  dearly 
seen.  ( 

The  mean  frequency  f  of  each  Doppler  spectrum  is 
calculated  by  the  equation, 

/-szyyp*  (6) 

where  Dk  is  the  power  of  the  spectrum  at  frequency  fk. 
The  partial  aliasing  of  each  of  all  the  8610  Doppler  spectra 
was  corrected  by  the  iteration  technique  until  /"converges 
to  /"*.  The  corresponding  spectra  in  Figs.  3  (a,b,c,d)  after 
aliasing  correction  were  replo'tted  in  Figs.  3  (a’,b',c',d'), 
where  each  spectrum  looks  much  more  symmetric  and 
much  closer  to  Gaussian  distribution  which  was  represent¬ 
ed  by  the  smooth  curve  in  the  same  figure.  The  compari¬ 
sons  between  the  mean  velocity  u  converted  from  the 
mean  frequency  f  of  each  spectrum  in  Figs.  3  (a,b,c,d) 
and  the  mean  velocity  u*obtained  from  the  corresponding 


Fig.  3.  (a,b,c,d)  Examples  of  power  spectra  obtained  by  finite- 
Fourier  transforms.  Partial  aliasings  are  clearly 
seen. 

(aib'.c'.d',)  Corresponding  power  spectra  of  fig.  3  fa.b.c.d) 
after  aliasing  correction  by  iteration  technique. 


spectrum  in  Figs.  3  (a'b’.c'd')  are  presented  in  Table  1, 
which  shows  that  partial  aliasing  of  Doppler  spectrum 
may  cause  the  velocity  measurement  error  as  high  as  80%. 

Partial  aliasing  is  due  to  the  aliasing  of  only  pan  of 
the  frequencies.  On  the  other  hand,  there  may  aiso  exist 

Table  1.  Comparison  of  Vertical  Velocity  v  Before  Aliasing  Cor¬ 


rection  and  C*  After  Aliasing  Correction 


Sample 

u  (cm/iec) 

u»  (cm/sec) 

error  1%) 

a 

15.4 

VS 

79 

b 

45.2 

62.4 

38 

c 

30.3 

31.6 

4.3 

d 

40.2 

57.3 

43 
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the  possibility  that  all  the  frequencies  of  the  spectrum  are 
aliased  by  the  same  number  of  cycles,  namely,  by  ruVA/ 
where  n  is  an  integer.  Also,  for  convenience  of  discussion, 
let’s  call  this  type  of  aliasing  as  total  aliasing  of  Doppler 
spectrum.  It  is  clear  that  in  the  case  of  total  aliasing,  the 
mean  frequency  J*  should  be  corrected  by  nNLf.  In 
principle,  the  vertical  velocity  can  be  any  one  of  the  fol¬ 
lowing  values, 

*(/*+nrVA/),  n  *  0,  ±1,  ±2, . . . 

Experience  from  conventional  measurements  suggests  that 
the  vertical  velocity  is  usually  in  the  order  smaller  than  1 
m/sec,  and  that  will  limit  n  to  be  either  0  or  ±1.  We  still 
have  to  tell  which  one  of  0,  +1,  and  -1  for  n  is  closest  to 
reality.  For  single  radar  beam  measurement,  the  possible 
solution  comes  from  the  comparison  between  the  present 
velocity  V*,  and  both  the  velocity  one  time  step  before, 
Vb*  and  the  velocity  one  time  step  after,  Vf.  If  both  V*  - 
V*  and  V*  -  V*  are  larger  than  \  A f),  then  V*  - 
^-.iVA/is  probably  closer  to  reality  than  V*.  On  the  con¬ 
trary,  if  both  Vp *  -  V*  and  V*  -  Vb *  are  smaller  than  -  $■ 
(^A /),  then  V*  +  ^ N&f  is  probably  better  than  V*. 
Otherwise,  Y*  is  probably  the  best  choice.  Or,  we  may 
also  compare  the  present  velocity  with  that  of  the  nearest 
upper  and  the  nearest  lower  height  ranges  and  do  the  same 
type  of  correction. 

For  three  or  more  spatially  separated  VHF  radar 
beam  measurement,  still  another  type  of  correction  of 
total  aliasing  is  possible.  The  following  section  is  an  ex¬ 
ample. 

Total  Aliasing  Correction  of  Vertical 
Velocities  Measured  by  Three  Spatially 
Separated  VHF  Radar  Beams 

The  data  evaluated  in  this  section  were  taken  with 
the  same  radar  by  the  same  institute  as  described  in  the 
previous  section,  starting  from  4PM,  October  28,  1981 
and  lasting  for  16  days.  The  radar  was  operated  either 
continuously  for  12  minutes  then  stopped  for  the  remain¬ 
ing  48  minutes  of  each  hour,  or  continuously  operating 
for  many  hours,  there  are  a  total  number  of  485  time 
intervals  of  operation  with  12  minutes  for  each,  and  the 
data  covers  1.35km-19km  height  range.  Three  antenna 
sets  are  arranged  as  shown  in  Fig.  4,  with  each  set  operat¬ 
ing  in  cyclic  order  for  1/3  seconds  then  stopping  for  2/3 
seconds  alternatively. 

The  raw  data  were  processed  by  Max-P!ancklnsii- 
tute  for  each  height  range  to  obtain  the  two  moments 
of  the  power  spectrum  of  radar  echo  signals  from  each 
radar  beam  of  each  12  minutes  time  period  by  the  auto¬ 
correlation  method. 

The  autocorrelation  function  of  a  coherently  detect- 
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Fig.  4.  Arrangement  of  3  spatially  separated  antenna  subsystems, 
where  ®  denotes  the  center  of  antenna  subsystem,  each 
containing  32  yagis,  and  the  average  transmitter  power  was 
■  20  lew. 

ed  radar  echo  signal  £(r)  is, 

C(r) *  (£(r)£*(r  +  r)> 

the  power  spectrum  of  the  signal  is, 

D (w)  *  yj/C(r)  exp  (-iur)dr 

here,  t  *  time,  cj  *  angular  frequency,  r  *  temporal  dis¬ 
placement  and  •  denotes  the  complex  conjugates.  The 
first  three  moments  m0,  mu  m7  of  the  spectrum  yield 
essential  parameters  of  the  radar  echo,  namely,  the  power 
P,  the  mean  frequency  J  and  the  spectral  width  a.  These 
are, 

P  =  mb  =  fD(oj)doj 
f  =  ^4  with  m,  =  fuD(u)d w 


0  -  with  m7  3  /wJD(w)cfc J. 

The  three  moments  can  also  be  deduced  directly  from  the 
autocorrelation  function  (Woodman  and  Guillen  [2] ) 
yielding 

P*C(,r0)  (r„*0) 

_ .  \I/(t  ) 

J  •  7  44-  (r.  *  first  time  lag) 

Ti 

f-jfr,)M(T8)}  )*  (7) 

r? 
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where  A{r)  is  the  amplitude  and  i£(r)  is  the  phase  of  the 
autocorrelation  function, 

C(r)  *  A{r)  exp  [i«4r(r)] 
or 

,  L  <2(r) 

*<r)',a”'  1A)-  w 

Since  the  arctangent  function  is  a  periodic  function 
with  period  2tr,  and  its  principal  value  is  defined  in  the 
range  between  -Hr  and  -rr,  so  the  phase  i It  given  by  Eq.  (8) 

can  be  aliased  from  its  true  value  by  2ntt,  n  =  1,2, _ 

This  is  just  the  case  of  total  aliasing,  its  possible  solution 
can  be  observed  through  comparisons  among  the  three 
phases  4/lt  <jj2,  i {/3  of  the  three  autocorrelation  functions 
of  the  three  respective  radar  echo  signals.  Since  the  se¬ 
paration  between  any  two  of  the  three  antenna  sets  is  less 
than  40  meters,  which  is  small  compared  with  the  scale 
of  the  atmospheric  motion  on  the  12  minutes  average,  it 
is  natural  to  expect  that  the  distributions  of  4/i  -  4>j  {i,j  = 
1,2.3)  should  follow  Gaussian  distribution  with  zero 
mean.  That  means  that  a  smaller  absolute  value  of  i ht  -  4>, 
is  always  more  probable  than  a  larger  absolute  value  of 
Wj  -  4sj-  This  principle  highlights  the  possibility  of  correct¬ 
ing  the  total  aliasing,  and  the  following  four  cases  as  de¬ 
monstrated  in  Fig.  5  ^a,b,c,d)  can  be  corrected  as  follows: 
Case  1.  In  Fig.  5a 

\4>,  -  4>k I  >  *,  1 4>r-  4/k\ >  *  and  -ir  <  J/k  <  -  Cor¬ 
rect  4/k  by  +2jt,  i.e.,  4>k  *  4/k  +  2jt,  then  both  l4/{~^k\  be¬ 
come  less  than  it.  We  can  see  that  is  probably  closer 
to  reality  than  4/k. 


Fig.  J.  Separations  of  three  phases  J(.  it/,  ik  {i.i.k  denote  1.2,3) 
of  the  three  autocorrelation  functions  of  three  radar  echo 

singnals. 


Case  2.  In  Fig.  5b 

With  similar  reason  of  Case  1,  we  may  correct  4/,  to 
3)  by  i b;-  2jt. 

Case  3.  In  Fig.  5c 

1 4>(  -  4>j\  >  it,  -  4/k\  >  tt  and  0  <  £,  <  y .  In  prin¬ 
ciple  we  may  correct  4/(  to  ^  by  » <j/t  -  2rr  to  make  both 
1 4>i  -  4>k  1  and  !yf-  4>j\  less  than  it.  However,  because  lv(.|  < 
y  and  lv,l  >  y.  1^1  seems  to  be  much  less  probable  than  i/i 
to  be  true.  So  instead  of  correcting  •!/.,  it  may  be  better 
to  correct  4>j  and  4>k  by  =  4ii  +  lit  and  1ik  =  +  2ir. 

Case  4.  In  Fig.  5d 

With  same  reason  of  Case  3,  we  may  correct  t b.  and 
4/k  by  =  4/j  -  2 it  and  vk  =  4>k-  2ir. 

As  examples,  the  means  and  variances  of  the  dis¬ 
tributions  of  the  phase  differences,  i/x  -  4t2,  4>\-  4>2  and 
4/ 2  -  4/ 2  °f  each  height  range  were  calculated  both  before 
and  after  aliasing  correction  and  the  results  were  present¬ 
ed  in  Figs.  6  -  11 .  Figs.  6-8  show  that  the  amplitude  of 
the  means  ■J/l  -  4>2  and  >//,  -  1 1/2  tend  to  be  too  large  for 
many  height  ranges  before  aliasing  corrections,  while  the 
amplitude  of  the  means  of  4/l  -  i^;,  ijr,  -  4/Jt  and  4>2  -  4>2 
are  generally  much  closer  to  zero  after  aliasing  correction, 
and,  Figs.  6  —  11  show  criteria  for  total  aliasing  correction 
do  make  sense.  Also,  Figs.  9-11  show  that  the  variances 
of  4>t  -  4>j  (/,/  =  1 ,2,3)  become  much  smaller  after  aliasing 
correction,  so  finer  distributions  for  4>f  were  obtained. 

It  is  important  to  make  one  point  clear;  there  is  no 
way  in  the  autocorrelation  analysis  to  detect  the  partial 
aliasing,  much  less  to  correct  it. 

Summary  and  Discussions 


Theoretically,  the  mean  frequency  /  of  a  power 
spectrum  as  defined  by  Eq.  (6)  and  that  deduced  directly 
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Height  range  (km) 

Figu  6.  j,  -  v,  at  each  height  range  before  aliasing  correcuon  (o) 
and  after  aliasing  correcuon  >•). 
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3.0  •  11.4  19.3 

Height  range  (km) 

Fig.  7.  -  v,  at  each  height  range  before  aliasing  correction  (o) 

and  after  aliasing  correction  (a). 

from  the  autocorrelation  function  as  given  by  Eq.  (7)  tend 
to  be  equal  to  each  other  as  the  time  step  A t  of  the  time 
series  approaching  to  zero,  and  the  total  time  interval  T 
approaching  to  infinity.  However,  both  A t  and  T  are 
finite  in  any  practical  case,  and  that  results  in  the  exist¬ 
ence  of  partial  aliasing  as  well  as  total  aliasing  of  frequen¬ 
cies  of  the  spectrum.  By  spectral  analysis  method  as  de¬ 
scribed  in  §2,  both  can  be  systematically  corrected  with 
confidence.  On  the  contrary,  there  is  no  way  by  autocor¬ 
relation  method  to  detect  the  partial  aliasing,  much  less 
to  correct  it.  Furthermore,  the  nonuaiform-amplification 
of  the  spectrum,  which  is  resulted  from  the  on-line  pre¬ 
integration  of  radar  echo  signals,  can  be  fixed  by  the  spec¬ 
tral  analysis  but  not  by  the  autocorrelation  method.  That 
underlines  the  shortcoming  of  the  autocorrelation  method 
proposed  by  Woodman  and  Guillen  [2]. 

The  iteration  method  to  correct  the  partial  aliasing 
problem  as  described  in  §2  is  very  powerful  to  eliminate 
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Height  range  (km) 

Fig.  3.  v,  -v,  at  each  height  range  before  aliaiing  correction  U) 
and  after  aliasing  correction  U). 
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Height  range  (km) 

Fig.  9.  Standard  deviation  (7  of  •*,  - 1 4,  distribution  at  each  height 
range  before  aliasing  correction  (o)  and  after  aliasing  cor¬ 
rection  (•). 

first  order  aliasing  of  spectra.  However,  one  should  be 
careful  that  it  might  not  be  applicable  for  second  and 
higher  order  aliasing. 

The  total  aliasing  correction  of  vertical  velocity 
measured  by  three  spatially  separated  VHF  radar  echoes 
was  based  on  our  belief  that  a  smaller  absolute  value  of 
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Height  range  (km) 

Fig.  10.  Standard  deviation  a  of  -  v,  distribution  at  each  height 
range  before  aliasing  correction  (3)  and  after  aliasing  cor¬ 
rection  (•). 


-  188  - 


■>  "J> ''  -  <  _•  >  _• 


V  _ 


I 

Aliasing  Problems  in  VHF  Radar  Signal 


1 _ i _ i - 1 - 1 - 1 — 

3.Q  11.4  19.8 

Height  range  (km) 

Fig.  11.  Standard  deviation  J  of  j/j  -  distribution  at  each  height 
range  before  aliasing  correction  (a)  and  after  aliasing  cor¬ 
rection  (a). 

•it.  -  i b.  is  always  more  probable  than  a  larger  value,  and 
we  think  that  this  also  holds  for  samples  which  may  be 
aliased  because  of  undersampling  a  time  series  varying 
because  of  a  mean  vertical  velocity  of  the  reflector.  As¬ 
sume  an  atmospheric  “reflector”,  which  has  an  extension 
larger  than  the  overlapping  Fresnel  zones,  or  antenna 
beams,  of  the  3  antennas.  The  reflector  shall  not  change 
its  shape,  but  shall  move  with  constant  vertical  speed.  We 
then  get  the  3  autocorrelation  functions  at  the  3  antennas 
equal  in  phase  and  amplitude.  The  correlation  time,  i.e., 
the  time  it  takes  for  the  amplitude  to  change,  is  very  long 
and  also  the  correlation  distance  is  very  large,  i.e.,  the 
signals  at  ail  3  antennas  are  similar.  Regardless  of  the  time 
when  the  samples  are  taken,  \pf  -  \j/.  =  0.  This  is  also  the 
case  for  any  vertical  velocity,  which  even  may  lead  to 
higher  order  aliasing. 

Now  assume  that  the  reflector  changes  its  shape  (say 
in  some  random  manner),  which  results  in  shorter  correla¬ 
tion  time  and  smaller  correlation  distance.  As  a  result,  the 
phases  at  tags  larger  than  the  correlation  time  are  changing 
randomly.  Also  the  phases  may  change  randomly  at  the  3 
antennas,  even  at  r  =  0,  because  the  spatial  correlation 
may  be  lost.  For  lags  shorter  than  the  correlation  time, 
the  phases  are  distributed  (say  Gaussian)  around  a  mean 
phase,  which  is  given  by  the  mean  vertical  velocity.  As¬ 
suming  that  the  mean  vertical  velocity  is  similar  overhead 
all  3  antennas,  the  mean  phases  must  be  similar.  The  dif¬ 


ferences  'h[  -  41/  are  distributed  around  zero.  The  width 
of  the  distribution  of  <hi  -  does  not  depend  on  the  mean 
vertical  velocity  rather  on  the  random  (turbulent)  velocity 
fluctuations.  Thus,  the  wider  this  distribution,  the  more 
variable  is  the  probed  atmospheric  region.  This  variability 
is  due  to  the  turbulent  fluctuations  of  the  region  of  the 
“reflector”  itself,  as  well  as  due  to  the  fact  that  this  re¬ 
flector  is  moving  horizontally. 

Then  we  propose  the  following  criterion  to  treat  the 
data.  Firstly,  we  assume  that.aiiasing  due  to  mean  vertical 
velocity  is  unlikely  (Nyquist  vertical  velocity  is  2.7  m/sec), 
but  it  may  occur  due  to  fluctuations.  There  then  is  a 
choice  to  select  among  good  data  which  will  not  be  cor¬ 
rected,  marginal  data  which  can  be  corrected,  and  bad 
data  which  have  to  be  corrected.  For  example,  it  can  be 
fair  that  only  those  sets  of  data  with  |  |  larger  than  a 

specified  variance  a  (say,  3  a  of  a  Gaussian  distribution) 
are  to  be  corrected  and  this  variance  a  should  be  predeter¬ 
mined  by  the  data  of  a  rather  quiet  time  period.  Of 
course.  Only  those  data  sets  which  indicate  a  sufficient 
signal-to-noise  ratio  shall  be  considered.  Still,  this  correc¬ 
tion  does  not  eliminate  any  aliasing  which  may  be  due  to 
a  large  vertical  velocity  (mean),  which  can  be  checked 
only  when  comparing  \[ilt  with  another  -  measured 

with  a  higher  sampling  rate. 

Additionally,  any  sets  of  •jti  -  \jj-  which  do  not  be¬ 
long  to  a  Gaussian  distribution,  e.g.,  rectangular  distribu¬ 
tion  in  case  of  strong  fluctuations,  should  be  taken  with 
great  care.  Instead  of  correcting,  a  more  acceptable  ap¬ 
proach  from  our  point  of  view  would  discard  those  data 
and  not  using  them  to  estimate  the  mean  vertical  velocity, 
but  still  use  the  criterion  to  determine  times  and  ranges 
which  were  strongly  variable  (i.e.  most  likely  turbulent). 

Acknowledgements 

One  of  the  authors  (F.S.  Kuo)  wants  to  thank 
EISCAT  and  Cornell  University  for  their  hospitality, 
where  part  of  this  work  is  carried  out,  and  to  thank  Na¬ 
tional  Science  Council  for  the  generous  support  of  his 
traveling  cost. 

References 

1.  Sato.  T.  and  R.F.  Woodman.  1982.  Spectral  parameter  estima¬ 
tion  of  CAT  radar  echoes  in  the  presence  of  fading  clutter. 
Radio  Science.  17,  =2,317-826. 

2.  Woodman.  R.F.  and  A.  Guillen.  1974.  Radar  observations  of 
winds  and  trubulence  in  the  stratosphere  and  mesosphere.  1. 
Atmos.  ScL.  31,  493-505. 


-  189  - 


F.S.  Kuo,  J.  Rattger,  J.K.  Chao,  C.H.  Liu  and  M.F.  Larsen 


VI @  ifi  £  'M  '< '&  M  !l  Jl  S  ibl  StJ  S  ¥  §iS  fi  S3  S 


IP  S  mm  •  J.  RSttger"  •  m  &  ft 


;ili  ?n  •  M.F.  Larsen* 


•  a2***$n»a* 

'*  EISCAT  Scientific  Association 


***** 


«k  « 

ias*«aa**a!E*«***sajM&as  •  xisutsazs  •  dittos 

ni55-*ffs#  •  SAifWiiaaawa^tittH*  •  •  «»it,«afe7»aw3‘ffr£itdi# 


Reprinted  from  Monthly  Weather  Review,  Vol.  Ill,  No.  10,  October  1983 

American  Meteoroiofical  Society 
Printed  m  U.  S.  A. 

Can  a  VHF  Doppler  Radar  Provide  Synoptic  Wind  Data? 

A  Comparison  of  30  Days  of  Radar  and  Radiosonde  Data 

Miguel  Folkmar  Larsen 

School  of  Electrical  Engineering,  Cornell  University,  Ithaca.  W  14853 
(Manuscript  received  4  January  1983,  in  final  form  16  June  1983) 

ABSTRACT 

A  number  of  experiments  have  shown  that  UHF  and  VHF  [Soppier  radars  can  make  “clear  air”  wind 
measurements  in  the  troposphere  and  lower  stratosphere,  even  in  the  presence  of  clouds  and  precipitation. 
Past  comparisons  of  radar  and  rawinsonde  profiles  for  a  single  day  have  shown  good  agreement  between  the 
two  within  the  limitations  imposed  by  the  spatial  separation  between  the  two  measurements.  However,  the 
accuracy  of  the  measurement  does  not  insure  that  the  data  are  synoptically  significant.  A  comparison  of  30 
days  of  radiosonde  geopotential  height  and  wind  data  and  radar  data  from  Alaska  was  carried  out  to  determine 
the  applicability  of  the  radar  data  to  synoptic  meteorology. 

Rawinsonde  and  radar  wind  time  series  at  six  heights  between  3.79  and  14.79  km  altitude  were  compared 
and  the  rms  differences  calculated.  The  two  independent  wind  measurements  were  also  compared  to  the 
geostrophic  wind  obtained  from  the  geopotential  height  fields  generated  by  applying  the  Cressman  objective 
analysis  scheme  to  data  from  five  radiosonde  nations  surrounding  the  radar  site.  The  two  independent  wind 
measurements  were  most  similar  with  a  difference  of  3-4  ms'1.  The  radar  and  balloon  winds  both  differed 
from  the  geostrophic  winds  by  5-6  m  s'1.  The  rms  differences  decreased  when  the  radar  winds  were  averaged 
over  longer  time  intervals.  The  cross  correlation  between  the  measured  and  geostrophic  winds  was  found  to 
be  —0.75  and  essentially  independent  of  height  over  the  altitude  range  studied. 


1.  Introduction  PROFS  (Prototype  Regional  Observing  and  Forecast- 

__  ing  System)  program  designed  to  improve  mesoscale 

The  number  of  sensitive  VHF  and  UHF  Doppler  forecasting  over  periods  of  a  day  or  less.  PROFS  in- 

radars  being  used  for  atmospheric  research  has  in-  eludes  a  network  of  three  VHF  radars  located  in  a 

creased  considerably  since  Woodman  and  Guillen  triangle  centered  at  Boulder,  Colorado.  Radar  wind 

(1974)  first  reported  on  their  measurements.  Radars  profiles  are  routinely  processed  and  transmitted  to  the 

operating  at  wavelengths  from  tens  of  centimeters  to  forecasting  office  once  an  hour  to  aid  in  short-range 

several  meters  measure  the  winds  from  the  Doppler  forecasting  (Hogg  et  al.,  1980;  Strauch,  1981;  Strauch 

shift  of  the  signal  backscattered  from  turbulent  van-  et  al.,  1982). 

ations  in  the  refractive  index.  Typical  height  and  time  The  relatively  low  cost  of  establishing  a  facility  is 
resolutions  are  of  the  order  of  hundreds  of  meters  and  part  of  the  appeal  of  the  radar  systems.  Also,  very  little 

several  minutes.  The  technique  and  some  of  the  early  maintenance  or  other  manual  intervention  is  required 
results  have  been  reviewed  by  Gage  and  Balsley  (1978),  once  the  facility  is  in  operation  (see,  e.g.,  Balsley  et 
Rottger  (1980)  and  Harper  and  Gordon  (1980).  al.,  1980).  The  present  study  was  undertaken  to  de- 

The  first  research  applications  of  the  long-wavelength  termine  if  a  VHF  radar  system  of  the  type  that  would 

radars  were  in  the  area  of  microscale  dynamics,  most  likely  be  used  operationally  can  provide  wind 

including  turbulence  in  the  troposphere  and  lower  data  that  is  synoptically  meaningful.  The  radar  provides 

stratosphere  and  wave  dynamics  at  scale  sizes  typical  a  very  accurate  wind  measurement,  i.e.,  a  very  accurate 

of  gravity  waves.  Over  the  past  few  years  there  has  measurement  of  the  velocity  of  the  turbulent  scatterers 

been  increasing  interest  in  using  the  radar  measure-  in  the  direction  parallel  to  the  beam.  The  attainable 

ments  for  synoptic  research  and  as  part  of  the  standard  accuracy  is  much  greater  than  what  would  be  useful 

observational  network.  Larsen  and  Rottger  (1982)  have  in  practice.  However,  the  measurement  may  be  con- 

reviewed  research  applications,  and  Balsley  and  Gage  taminated  with  high-frequency  subsynoptic-scale  noise, 

(1982)  have  discussed  some  of  the  considerations  in-  making  the  data  unusable.  Also,  the  VHF  radars  typ- 

volved  in  establishing  an  operational  observing  system,  ically  operate  with  large  fixed  dipole  arrays  that  look 

Carlson  and  Sundararaman  (1982)  have  made  strong  within  a  few  degrees  of  vertical.  Thus  scans  of  large 

arguments  for  the  economic  feasibility  of  implementing  horizontal  areas  to  determine  the  wind  representative 

such  a  network.  One  of  the  most  ambitious  operational  of  some  area  greater  than  a  predetermined  horizontal 

projects  utilizing  the  radar  technique  is  part  of  the  scale  size  are  not  possible.  However,  it  may  be  that 
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averaging  high  time-resolution  wind  profiles  can  pro¬ 
vide  a  suitable  filtering.  The  next  sections  will  inves¬ 
tigate  the  amount  of  averaging  that  is  necessary. 

A  comparison  limited  to  the  rawinsonde  and  radar 
winds  would  not  reveal  if  the  differences  are  due  to 
the  peculiarities  of  one  measurement  system  or  the 
other.  Therefore,  the  two  independent  wind  measure¬ 
ments  were  also  compared  to  the  geostrophic  wind 
calculated  from  a  grid  of  geopotential  heights  derived 
from  nearby  radiosonde  measurements.  The  geo¬ 
strophic  wind  then  provides  an  independent  standard 
for  comparison. 

2.  Synoptic  representativeness 

The  general  meaning  of  synoptic  representativeness 
is  dear,  based  on  the  discussion  above,  but  a  specific 
quantitative  comparison  was  needed  for  this  study. 
The  problem  is,  in  fact,  closely  related  to  the  problems 
of  data  assimilation  and  objective  analysis  for  nu¬ 
merical  models  (see  e.g.,  Bengtsson  et  al..  1981).  The 
latter  attempts  to  treat  a  set  of  observations  in  such  a 
way  that  high-frequency  waves  that  lead  to  numerical 
instabilities  are  eliminated.  A  balance  condition  of  one 
kind  or  another  is  used  to  accomplish  this,  the  simplest 
being  a  geostrophic  balance  (e.g.,  Cressman,  1959)  and 
the  most  involved  being  the  balance  derived  from  the 
model  dynamics  and  termed  normal  mode  initializa¬ 
tion  (e.g.,  Daley,  1981). 

In  this  study  I  will  compare  radar  wind  data  taken 
over  a  30-day  period  from  1  March  to  1  April  1979 
and  radiosonde  data  taken  at  the  standard  observa¬ 
tional  hours  during  the  same  period.  The  radar  is  lo¬ 
cated  at  Poker  Flat,  Alaska  [see  Balsley  el  al.  (1980) 
for  a  description],  and  the  nearest  radiosonde  station 
is  at  Fairbanks,  35  km  southeast  of  the  radar.  The 
radar  measures  one  complete  wind  profile  every  4  min 
between  the  surface  and  14.79  km,  but  these  high  time- 
resolution  measurements  are  likely  to  be  contaminated 
by  high-frequency  meteorological  '‘noise."  In  order  to 
test  whether  or  not  the  radar  can  provide  useful  syn¬ 
optic-scale  wind  information,  comparisons  were  made 
between  the  geostrophic  wind  derived  by  applying  an 
objective  analysis  scheme  to  the  geopotential  height 
data  from  five  surrounding  radiosonde  stations  and 
the  radar  winds  averaged  over  various  time  intervals. 

A  test  of  geostrophy  may  not  be  the  best  test  of 
synoptic  representativeness,  but  a  simple,  yet  useful, 
comparison  was  needed  for  the  purpose  of  this  initial 
study.  Also,  the  differences  between  the  radar,  rawin¬ 
sonde  and  geostrophic  winds  are  of  interest  from  the 
point  of  view  of  the  multivariate  objective  analysis 
schemes  that  use  geostrophy  to  relate  height-height 
autocorrelations  to  the  height-wind  cross  correlations 
(e.g.,  Rutherford,  1972:  Schlatter,  1975;  Schlatter  et 
al.  1976:  Bergman,  1979;  Bengtsson  et  al.,  1981:  Lor- 
enc.  1981).  Generally,  the  wind  measurements  are  used 
to  provide  information  on  the  geopotential  height  gra¬ 
dients  and  these  are  used  to  formulate  the  initial  state 


for  a  model  integration.  The  results  of  this  study  show 
that  the  radar  and  rawinsonde  measurements  are  most 
similar  in  the  rms  (root-mean-square)  sense,  and  either 
wind  measurement  differs  from  the  geostrophic  wind 
by  a  larger  amount  than  the  difference  between  the 
two  wind  measurements. 

I  will  describe  the  radar  and  radiosonde  data  sets 
in  the  next  section.  Radar  and  rawinsonde  wind  mea¬ 
surements  are  compared  in  the  following  section.  Next, 
radar  and  rawinsonde  measurements  are  compared  to 
the  geostrophic  winds  derived  from  the  geopotential 
height  data.  The  last  section  contains  a  general  dis¬ 
cussion  of  the  implications  of  the  results. 

3.  Description  of  the  data  sets 

During  the  spring  of  1979,  the  Poker  Rat  MST  radar 
operated  with  a  100  m  X  100  m  fixed  dipole  array. 
The  antenna  size  has  since  been  expanded  to  200  m 
X  200  m.  The  horizontal  wind  vector  was  measured 
using  two  beams  pointed  15°  off  vertical,  one  with  an 
azimuth  pointing  close  to  north  and  one  pointing  close 
to  east.  Corrections  were  applied  to  the  data  to  get  the 
exact  zonal  and  meridional  wind  components.  One 
complete  wind  profile  was  obtained  every  4  min  with 
a  height  resolution  of  2.2  km.  In  the  troposphere  and 
lower  stratosphere  wind  data  were  routinely  available 
at  six  consecutive  heights  beginning  at  3.79  km  altitude. 
All  altitudes  for  radar  data  referred  to  from  here  on 
are  at  the  center  of  the  range  gate.  The  wind  mea¬ 
surement  is,  in  some  sense,  an  average  over  a  2.2  km 
height  range  weighted  by  the  intensity  of  the  turbulent 
layers  in  the  volume  illuminated  by  the  radar  beam, 
and.  as  Sato  and  Fukao  (1982)  have  pointed  out.  having 
a  number  of  turbulent  layers  withn  the  range  gate  can 
lead  to  errors  in  the  measurement  if  there  is  a  strong 
shear  across  the  range  gate.  However,  such  an  effect 
will  not  be  considered  further  here. 

The  signal  voltages  measured  by  the  receiver  are 
Fourier  transformed  to  obtain  the  Doppler  spectrum. 
The  frequency  spectrum  measured  by  scattering  from 
a  turbulent  process  will  be  Gaussian  for  a  normal  ran¬ 
dom  process  as  discussed  by  Woodman  and  Guillen 
(1974)  and  others.  The  accuracy  with  which  the  velocity 
of  the  turbulent  scatterers  in  the  radar  volume  can  be 
determined  is  then  limited  only  by  the  rate  at  which 
the  signal  is  sampled  and  the  accuracy  with  which  the 
Gaussian  can  be  fit  to  the  spectrum  (see  e.g..  Fukao 
et  al.,  1982).  Accuracies  better  than  1  eras'1  can  be 
achieved  in  practice,  but  there  is  usually  no  need  to 
measure  the  horizontal  wind  components  that  accu¬ 
rately.  The  error  for  the  radar  data  discussed  here  is 
~  1  m  s'1-.  When  all  three  wind  components  are  mea¬ 
sured  simultaneously,  the  signals  from  the  vertical 
beam  are  coherently  integrated  for  a  longer  time  than 
the  signals  from  the  off-vertical  beams  to  attain  the 
much  higher  accuracy  needed  for  vertical  velocity 
measurements. 
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The  horizontal  vector  winds  were  derived  from  the 
off-vertical  wind  components  with  the  assumption  that 
the  vertical  velocity  was  negligible.  In  most  cases  such 
an  assumption  is  justified,  as  more  recent  data  sets 
have  shown  when  all  three  wind  components  were 
measured  (T.  Riddle,  Aeronomy  Laboratory,  NOAA, 
private  communication,  1982).  However,  there  are 
times  when  the  vertical  velocity  can  make  a  significant 
contribution  over  a  period  spanning  several  consecutive 
profiles.  In  order  to  minimize  any  such  effect,  the  radar 
wind  data  were  averaged  over  an  hour.  The  time  series 
of  the  hourly-average  zonal  and  meridional  wind  com¬ 
ponents  are  plotted  in  Fig.  1  taken  from  Larsen  et  al. 
(1982). 

The  radiosonde  data  provided  by  the  National  Cli¬ 
matic  Center  are  from  the  five  nearest  National 
Weather  Service  (NWS)  radiosonde  stations  as  shown 
in  Fig.  2.  Only  balloon  launches  at  0000  and  1200 
GMT  are  available  for  the  period  discussed  here.  No 


Fio.  2.  Map  showing  the  location  of  the  five  radiosonde  stations 
closest  to  the  site  of  the  Poker  Flat  radar. 
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Fig.  1.  Original  data  set  of  40  days  of  hourly  values  measured  by 
the  Poker  Rat.  Alaska  MST  radar  beginning  on  23  February  and 
ending  on  5  April  1979.  The  top  half  of  the  figure  is  for  the  zonal 
wind  and  the  bottom  half  for  the  meridional  wind  component.  Only 
the  period  from  1  March  to  1  Apnl  was  used  in  the  study. 


supplemental  launches  are  included  in  the  data  set. 
The  radiosonde  station  closest  to  Poker  Flat  is  at  Fair¬ 
banks,  35  km  southeast  of  the  radar  site.  Anchorage, 
McGrath,  Barter  Island  and  Yakutat  are  455, 476,  572 
asd  764  km  from  the  radar  site,  respectively. 

The  balloon  wind  data  were  interpolated  linearly 
between  the  two  nearest  significant  levels  to  the  altitude 
corresponding  to  the  center  of  each  of  the  radar  range 
gates.  The  interpolation  was  applied  independently  to 
the  wind  speed  and  the  wind  direction.  The  pressure 
at  the  seven  different  radar  heights  was  determined 
from  the  Fairbanks  radiosonde  data.  The  formula  for 
a  hydrostatic  atmosphere  with  a  constant  lapse  rate 
was  then  used  to  interpolate  the  pressure  between  the 
two  nearest  significant  levels,  and  the  same  formula 
was  applied  to  the  radiosonde  data  from  the  other  four 
stations  to  find  the  height  of  the  pressure  surface  in¬ 
tersecting  the  center  of  each  of  the  radar  range  gates. 

4.  Comparison  of  radar  and  rawinsonde  wind  mea¬ 
surements 

The  rawinsonde  and  radar  winds  for  the  period  from 
1  March  to  1  April  1979  are  plotted  in  Fig.  3  for  seven 
consecutive  heights.  The  heavy  solid  line  is  the  rawin¬ 
sonde  measurement  and  the  light  dashed  line  the  radar 
measurement.  The  signal-to-noise  ratio  at  16.99  km 
is  too  low  to  derive  a  meaningful  signal  from  the  echoes, 
but  the  time  series  at  this  height  have  been  plotted  for 
reference,  nonetheless.  The  radar  winds  are  1  h  av¬ 
erages,  as  mentioned  previously,  while  the  rawinsonde 
data  are  from  a  single  ascent.  The  rms  difference  be¬ 
tween  the  two  measurements  is  indicated  by  the  quan¬ 
tity  labeled  SIG  (sigma). 

The  two  curves  follow  each  other  very  closely  until 
day  78  but  deviate  significantly  from  each  other  after 
that.  The  rms  differences  range  from  a  minimum  of 
3.21  ms'1  to  a  maximum  of  9.30  ms'1  in  the  zonal 
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Fig.  3.  Comparison  of  1  h  average  radar  winds  (light  dashed  line) 
and  rawinsonde  wind  measurements  (heavy  solid  line)  at  seven  suc¬ 
cessive  heights  during  March  1979.  The  rms  (root-mean-square)  dif¬ 
ferences  are  shown  at  bottom  right-hand  side  of  each  bos. 


the  difference  is  due  to  small-scale  temporal  and  spatial 
variability  and  if  the  winds  are  statistically  well  be¬ 
haved.  Indeed,  that  is  the  case.  Fig.  5  shows  the  com¬ 
parison  between  the  24  h  average  rawinsonde  and  radar 
winds.  The  radar  average  consists  of  24  h  of  wind 
profiles  taken  every  4  min.  The  rawinsonde  average 
consists  of  three  wind  profiles,  one  taken  at  either  0000 
or  1 200  GMT  and  the  profile  preceding  and  following 
it.  The  rms  differences  have  decreased  to  a  maximum 
of  6.52  m  s’1  in  the  zonal  component  and  a  maximum 
of  7.33  m  s'1  in  the  meridional  component. 

A  comparison  of  the  12  h  average  radar  winds  and 
the  winds  measured  by  individual  rawinsonde  ascents 
was  also  made.  In  other  words,  the  averaging  interval 
was  increased  for  the  radar  when  compared  to  Fig.  3 
but  not  for  the  rawinsonde.  The  rms  differences  still 
decreased  significantly  when  compared  to  the  case  of 
the  1  h  averages.  The  difference  between  a  1  h  and  a 
12  h  average  of  the  radar  winds  was  primarily  that 
many  of  the  large  spikes  in  the  radar  data  have  been 
eliminated.  The  original  time  series  shown  in  Fig.  I 
has  quite  a  few  large  spikes.  Some  of  them  are  isolated, 
but  a  number  of  them  occur  during  periods  of  high 
variability,  indicating  that  they  are  likely  to  be  real. 
Averaging  decreases  the  large  errors  that  arise  due  to 
this  effect.  The  rawinsonde  winds  do  not  exhibit  any 
of  these  spikes  in  spite  of  the  fact  that  they  are  not 
averaged  over  long  time  intervals  as  the  radar  winds 
are.  However,  the  profiles  are  routinely  edited  to  delete 
wind  data  at  any  heights  characterized  by  large  changes 
in  wind  direction  or  wind  speed  over  1-2  min  intervals 
during  the  ascent.  An  attempt  was  made  to  subject 
the  radar  data  to  the  same  type  of  screening.  The  exact 


component  and  from  a  minimum  of  4 .66  ms-1  to  a 
maximum  of  12.73  m  s_l  in  the  meridional  compo¬ 
nent.  The  variation  of  sigma  with  height  follows  the 
variation  of  the  mean  wind  for  the  period  as  a  function 
of  height.  The  mean  wind  has  been  plotted  in  Fig.  4. 
Schlatter  (1975)  and  Jasperson  (1982)  also  found  that 
the  difference  between  the  objectively  analyzed  wind 
and  the  measured  wind  varied  as  the  mean  wind  profile. 

The  difference  in  the  winds  measured  at  the  two 
locations  with  the  two  different  techniques  can  be  at¬ 
tributed  to  a  combination  of  three  factors.  First,  there 
are  errors  in  the  wind  measurements;  second,  there  is 
variability  in  the  wind  over  the  spatial  separation  be¬ 
tween  the  two  sites;  and  third,  there  is  temporal  vari¬ 
ability.  The  rawinsonde  takes  approximately  1  h  to 
ascend  through  the  height  range  where  reliable  mea¬ 
surements  can  be  made  by  the  radar.  The  radar  mea¬ 
surements  presented  here  are  1  h  averages,  but  the 
rawinsonde  will  spend  only  a  few  minutes  within  any 
one  of  the  radars  range  gates. 

It  should  be  possible  to  decrease  the  rms  differences 
between  the  two  wind  measurements  by  averaging  the 
measurements  over  longer  time  intervals  if  some  of 
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Fig.  4.  Mean  wind  profile  for  the  32-dav  period  based  on  the 
radar  wind  measurements.  The  zonal  and  meridional  wind  com¬ 
ponents  were  averaged  separately  to  determine  the  curves. 
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Fig.  5.  Comparison  of  24  h  average  radar  winds  (light  dashed 
line)  and  averages  of  three  successive  rawinsonde  wind  measurements 
(heavy  solid  line). 


criteria  could  not  be  applied  since  1-2  min  of  balloon 
ascent  correspond  to  approximately  300-600  m  of  al¬ 
titude.  The  height  resolution  of  the  radar  measurements 
is  2.2  km.  Thus,  the  rawinsonde  criteria  were  linearly 
extrapolated  by  multiplying  the  change  over  600  m 
by  the  appropriate  factor  to  get  the  corresponding 
change  over  2.2  km.  There  is  little  doubt  that  the  two 
procedures  are  not  equivalent.  In  any  case,  applying 
the  error  criteria  did  not  change  the  results  significantly. 

Figs.  3  and  5  both  show  a  similar  behavior  in  one 
respect.  The  agreement  between  the  radar  and  rawin¬ 
sonde  wind  measurements  is  quite  good  until  day  78. 
After  that  a  large  discrepancy  arises.  One  would  suspect 
that  most  of  the  rms  difference  between  the  two  mea¬ 
surements  is  attributable  to  the  period  from  day  78 
until  the  end.  I  will  show  later  than  the  statistics  are 
quite  different  for  the  first  half  and  the  last  half  of  the 
month  and  discuss  possible  reasons  for  such  an  effect. 

5.  Comparison  of  measured  winds  and  geostrophic 
winds 

The  Cressman  analysis  is  an  objective  analysis 
scheme  used  to  interpolate  data  from  irregularly  spaced 
grid  points  to  a  regular  grid  suitable  for  analysis  or 


input  to  a  numerical  model.  Given  a  series  of  obser¬ 
vations  /'  where  the  index  i  ranges  over  the  number 
of  observations  that  influence  the  grid  point,  the  value 
of  /at  the  grid  point  ft  will  be  given  by 

/i-/J+  2^/',.  (1) 

f-l 


The  quantity  f%  is  an  initial  guess  of  the  value  of / at 
the  grid  point.  The  initial  value  will  be  modified  by 
the  observed  changes  from  the  initial  guesses  at  the 
observation  points  given  by  f,  =  f°  -/?.  The  actual 
observed  value  is  /?.  The  weighting  function  is  given 
by  C{  =  n~'.  (R2  -  ^/(R2  +  D2),  where  R  is  the 
radius  of  influence  and  D  is  the  distance  from  the 
observation  point  to  the  grid  point.  The  weight  is  zero 
if  D  is  greater  than  R.  The  distance  D  was  calculated 
from  the  formula  given  by  Schlatter  (1975) 

D  -  rE j(0,  -  cos2^(9?  +  0,)J 


1/2 


(2) 


where  rE  is  the  radius  of  the  earth,  8  is  latitude  and  <t> 
is  longitude.  The  form  of  the  objective  analysis  scheme 
used  here  has  been  described  in  detail  by  Kruger  (1969). 

The  geostrophic  wind  components  were  derived 
from  the  observed  values  of  the  geopotential  height  at 
the  five  radiosonde  stations  shown  in  Fig.  2.  Values 
were  interpolated  to  a  grid  with  a  spacing  of  100  km 
and  aligned  along  the  north-south  and  east-west  axes. 
The  central  grid  point  was  positioned  to  coincide  with 
the  location  of  the  Poker  Flat  radar.  Typically,  the 
interpolation  procedure  would  be  confined  to  fixed 
pressure  surfaces.  However,  the  radar  measures  the 
velocity  at  a  fixed  height  irrespective  of  pressure.  Thus, 
the  initial  guesses  of  the  geopotential  height  at  the 
observation  points  had  to  be  derived  in  a  slightly  dif¬ 
ferent  way.  The  pressure  at  the  height  of  the  seven 
different  radar  ranges  was  calculated  based  on  the  Fair¬ 
banks  sounding.  The  height  of  the  corresponding  pres¬ 
sure  surface  at  the  four  surrounding  radiosonde  stations 
was  calculated,  and  the  derived  height  values  at  a  given 
station  were  then  averaged  for  the  entire  30-day  period 
during  which  data  were  available.  Therefore,  the  initial 
or  climatological  value  of  the  height  at  a  given  radio¬ 
sonde  station  is  the  average  geopotential  height,  for 
the  period,  of  the  average  pressure  surface  correspond¬ 
ing  to  the  height  of  the  radar  range  gate  at  Poker  Flat. 

The  24  h  average  radar  winds  are  plotted  in  Fig.  6 
together  with  the  geostrophic  wind  calculated  from  the 
height  field  obtained  by  applying  the  procedure  de¬ 
scribed  above.  The  rms  differences  are  larger  than  those 
between  the  radar  and  rawinsonde  wind  measurements. 
Also,  the  errors  are  much  larger  in  the  meridional 
component  than  in  the  zonal  component.  It  appears 
that  most  of  the  contribution  to  the  large  rms  differ¬ 
ences  come  from  the  second  half  of  the  time  series. 
Fig.  7  is  similar  to  Fig.  6  but  shows  the  24  h  average 
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Flo.  6.  Comparison  of  24  h  average  radar  winds  (light  dashed 
line)  and  the  geostrophic  wind  (heavy  solid  line).  The  geostrophic 
wind  was  calculated  from  a  grid  of  geopotential  height  values  derived 
using  the  Cressman  objective  analysis  scheme  with  a  radius  of  in¬ 
fluence  of  7?0  km. 

rawinsonde  winds  as  opposed  to  the  24  h  average  radar 
winds.  The  rms  differences  between  each  of  the  wind 
measurements  and  the  geostrophic  winds  are  com¬ 
parable. 

The  behavior  of  the  differences  between  the  radar 
winds  and  the  geostrophic  winds  is  summarized  in  Fig. 
8a  and  8b.  The  error  decreases  as  the  averaging  interval 
is  increased  in  both  the  zonal  and  meridional  com¬ 
ponent.  Although  the  error  decreases  significantly  when 
the  averaging  interval  is  increased  from  I  hour  to  12 
hours,  particularly  in  the  meridional  component,  the 
reduction  is  relatively  small  when  the  averaging  interval 
is  increased  further.  The  error  is  much  larger  in  the 
meridional  component  than  in  the  zona!  component. 
The  overall  shape  of  the  curves  is  very  similar  to  the 
shape  of  the  mean  wind  profile  shown  in  Fig.  4. 

The  rms  differences  between  the  three  24  h  averages 
are  shown  in  Fig.  9a  for  the  zonal  component  and  9b 
for  the  meridional  component.  There  is  very  little  dif¬ 
ference  between  the  three  for  the  zonal  component, 
although  the  radar  and  balloon  wind  measurements 
are  most  similar  at  the  two  highest  altitudes.  The  radar 
and  balloon  winds  are  most  similar  for  the  meridional 
component  except  at  the  two  lowest  altitudes.  It  appears 


that  in  general  the  differences  between  the  two  inde¬ 
pendent  wind  measurements  are  smaller  than  the  dif¬ 
ferences  between  the  measured  and  geostrophic  winds, 
but  the  changes  are  so  small  that  it  is  difficult  to  justify 
any  inferences  on  the  basis  of  these  results.  However, 
when  the  first  and  second  half  of  the  data  sets  are 
treated  separately,  the  relationship  is  much  clearer. 

6.  Comparison  of  first  and  second  half  of  data  set 

Close  examination  of  Figs.  3.  5,  6  and  7  shows  that 
most  of  the  contribution  to  the  rms  differences  between 
the  radar  and  rawinsonde  winds,  the  radar  winds  and 
the  geostrophic  wind,  and  the  rawinsonde  winds  and 
the  geostrophic  winds  derives  from  the  time  period 
between  day  78  and  day  92,  the  end  of  the  period 
when  data  are  available.  The  time  series  were  divided 
into  two  segments,  one  from  I  March  to  16  March 
and  the  other  from  17  March  to  1  April.  The  rms 
differences  were  calculated  for  each  half  for  1  and  24 
h  averages.  The  large  deviations  occur  later  than  17 
March,  but  the  data  set  was  divided  in  half  so  as  not 
to  change  the  statistics  simply  due  to  differences  in 
the  number  of  samples  in  each  set. 

Table  1  shows  that  all  errors  are  comparable  in  the 
second  half  of  the  data  set  for  both  the  1  and  24  h 
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Fig.  7.  As  in  Fig.  6  but  for  an  average  of  three  successive  ra- 
ninsonde  measurements  instead  of  the  radar  wind  measure  me  rj 
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Zone!  component 
Radar  wind*  v*  Cressman 
analysis  R*  750 


RMS  difference  (m/s) 

Fig.  8.  (a)  Summary  plot  of  the  rms  differences  between  the  zonal 
component  of  the  radar  winds  and  geostrophic  winds  for  various 
averaging  intervals. 

averages.  However,  the  relationships  that  emerged 
when  the  time  period  was  treated  as  a  whole  are  much 
dearer  in  the  first  half.  When  a  1  h  average  of  the 
radar  data  is  compared  to  a  single  rawinsonde  ascent, 
the  differences  are  3. 19  m  s*1  in  the  zonal  component 
and  4.03  ms*1  in  the  meridional  component.  The 
differences  between  the  radar  and  geostrophic  winds, 
and  the  rawinsonde  and  geostrophic  winds  are  nearly 
identical  and  1. 5-2.0  m  s-1  greater  than  the  difference 


Meridional  component 
Radar  winds  vs  Cressman 
analysis  R*  750 


RMS  difference  (m/s) 


RMS  differences  for 
24  hr  averages 
Zonal  component 


RMS  difference  (m/ s) 

Fig.  9.  (a)  Summary  plot  of  the  rms  differences  for  the  radar/ 
geostrophic  wind,  radar/rawinsonde  and  rawinsonde/geostrophic 
wind  comparisons  for  the  zonal  component. 

between  the  radar  and  rawinsonde  winds.  Averaging 
for  24  h  decreases  the  difference  between  the  two  in¬ 
dependent  wind  measurements  to  2.21  and  3.10  m 
s*'.  The  difference  between  the  measured  winds  and 
the  geostrophic  winds  is  still  1 .5-2.0  m  s*1  greater  than 
the  difference  between  the  two  wind  measurements. 

Fukao  e:  at.  ( 1 982)  compared  data  taken  with  the 
Arecibo  430  MHz  radar  and  rawinsonde  data  from 
San  Juan,  Puerto  Rico.  The  separation  between  the 
two  sites  is  ~75  km.  Their  radar  wind  profiles  were 
averages  over  30  min.  but  their  data  taking  scheme 
did  not  involve  long-time  sequences  such  as  those  in 
the  Poker  Flat  data.  Consequently,  they  could  not  av- 


RMS  dif ferences  ‘or 
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Meridional  comoonenf 
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Fig.  S.  (bl  As  in  Fig.  8a  but  for  the  meridional  wind  component. 


Fig.  9.  (bl  As  in  Fig.  9a  but  for  the  meridional  wind  component. 
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Table  1.  Root  mean  square  differences. 


3/1/79-  3/17/79- 

3/16/79  4/1/79 

U  V  U  V 


I -hour  averages 


Radar  vs.  Rawinsonde 

3.19 

4.03 

8.45 

10.70 

Radar  vs.  Cressman  Analysis 
Rawinsonde  vs.  Cressman 

5.67 

5.57 

8.40 

12.37 

Analysis 

6.54 

5.30 

6.08 

8.84 

24-hour  averages 

Radar  vs.  Rawinsonde 

2.21 

3.10 

5.59 

6.44 

Radar  vs.  Cressman  Analysis 
Rawinsonde  vs.  Cressman 

4.66 

4.60 

5.79 

8.95 

Analysis 

5.29 

4.60 

5.20 

8.22 

erage  over  longer  time  intervals.  A  total  of  26  separate 
rawinsonde  ascents  were  used  in  the  comparison.  They 
found  a  distinct  variation  in  the  magnitude  of  the  rms 
differences  above  and  below  the  tropopause.  The  tro¬ 
pospheric  value  was  4.9  m  s"1  and  the  stratospheric 
value  was  3.3  m  s'1.  The  difference  that  they  calculated 
for  the  stratosphere  is  in  good  agreement  with  the  value 
found  in  this  study  when  the  1  h  averages  of  the  radar 
winds  were  compared  to  single  rawinsonde  ascents. 
However,  there  was  no  indication  in  this  study  of  a 
change  in  the  magnitude  of  the  differences  when  cross¬ 
ing  from  the  troposphere  into  the  stratosphere,  al¬ 
though  there  was  a  variation  that  followed  the  mean 
wind  profile. 

Fukao  et  al.  (1982)  estimated  the  error  of  the  ra¬ 
winsonde  measurements  to  be  1 .9  m  s'1  for  the  average 
conditions  during  their  experiment  although,  on  oc¬ 
casion,  the  error  was  as  large  as  7.1  m  s_l  when  the 
wind  speeds  were  large.  They  attributed  the  differences 
in  the  stratosphere  to  errors  in  the  rawinsonde  mea¬ 
surements  and  the  larger  differences  in  the  troposphere 
to  a  combination  of  the  rawinsonde  errors  and  in¬ 
creased  variability  below  the  tropopause.  The  difference 
between  the  two  wind  measurements  in  this  study  was 
found  to  be  3-4  ms-1.  However,  it  cannot  be  attributed 
to  the  error  in  the  rawinsonde  measurement  since  the 
difference  between  both  wind  measurements  and  the 
geostrophic  wind  was  about  the  same,  indicating  that 
the  difference  is  due  to  the  spatial  and  temporal  vari¬ 
ability  at  scales  comparable  to  the  separation  between 
the  two  sites  and/or  the  averaging  interval. 

Jasperson  ( 1982)  studied  mesoscale  wind  variability 
using  a  series  of  balloons  tracked  by  the  METRAC 
positioning  system.  He  was  able  to  study  the  variability 
at  the  time  lags  from  30  min  to  5  h  and  spatial  sep¬ 
arations  from  20  m  to  20  km.  Separations  as  large  as 
35  km  were  not  studied,  but  Jasperson’s  results  can 
be  extrapolated  to  estimate  the  variability  correspond¬ 
ing  to  the  separation  between  Fairbanks  and  Poker 
Rat.  The  derived  value  is  close  to  3  m  s'1  when  a 
power  law  of  the  form  given  by  his  Eq.  (1)  is  used 
together  with  the  values  that  he  measured  at  4.415 


and  20.910  km  separation.  That  is  in  good  agreement 
with  the  differences  between  the  I  h  average  radar 
winds  and  the  rawinsonde  winds  in  this  study.  Jas¬ 
person’s  (1982)  results  also  suggest  that  the  differences 
in  the  rms  values  characteristic  of  the  troposphere  and 
stratosphere  found  by  Fukao  et  al.  (1982)  could  be 
due  to  the  effect  of  the  changing  separation  between 
the  balloons  and  the  radar  observation  site  above  and 
below  the  tropopause.  Their  Fig.  I  shows  the  balloon 
trajectories  for  the  26  individual  ascents.  Due  to  the 
prevailing  stratospheric  easterlies  during  August  and 
September  at  18°N,  the  trajectories  are  such  that  the 
rawinsondes  begin  to  approach  the  radar  site  once  they 
get  into  the  stratosphere.  The  effect  of  decreasing  vari¬ 
ability  with  decreasing  spatial  separation  would  be  to 
decrease  the  difference  between  the  two  measurements. 

7.  Effect  of  varying  the  radius  of  influence 

All  of  the  curves  so  far  have  shown  results  for  a 
Cressman  objective  analysis  scheme  with  a  radius  of 
influence  of  750  km.  This  particular  radius  was  chosen 
because  it  was  the  optimum  value  when  the  entire 
height  range  where  radar  data  were  available  was  con¬ 
sidered.  Figs.  10a  and  1 0b  show  the  rms  differences 
between  the  24  h  average  radar  winds  and  the  geo¬ 
strophic  winds  as  a  function  of  the  radius  of  influence. 
The  750  km  value  gives  the  best  result  over  all,  although 
at  the  two  lowest  altitudes  the  best  result  is  achieved 
when  R  =  1000  km.  However,  the  difference  in  the 
errors  when  the  larger  value  is  used  is  very  small  below 
the  tropopause  but  can  amount  to  1-2  m  s'1  at  the 
tropopause  and  in  the  lower  stratosphere.  A  radius  of 
influence  of  750  km  also  gives  the  best  result  for  the 
meridional  component,  but  the  error  is  much  less  de¬ 
pendent  on  the  radius  chosen  in  this  component  than 
in  the  zonal  component. 

The  Gandin  objective  analysis  scheme  described  by 
Kruger  ( 1 969)  was  also  tested.  The  basic  difference  is 
in  the  weighting  function.  While  the  Cressman  analy¬ 
sis  uses  a  parabolic  weighting  function,  the  Gandin 
method  uses  a  Gaussian  weighting,  and  both  schemes 
are  univariate.  1  will  not  present  the  results  in  graphic 
form  since  the  difference  in  all  cases  was  only  a  fraction 
of  a  meter  per  second,  less  than  the  accuracy  of  the 
radar  measurements  during  the  spring  of  1979. 

The  Cressman  analysis  is  often  applied  in  a  series 
of  successive  scans  (Cressman,  1959:  Otto-Bliesner  et 
al.,  1977).  A  large  radius  of  influence  is  chosen  initially 
so  as  to  include  information  from  many  stations  sur¬ 
rounding  the  grid  point.  The  radius  is  then  decreased 
successively,  until  finally  only  the  nearest  stations  affect 
the  analyzed  value.  The  successive  scans  method  was 
also  applied  to  the  Alaska  radiosonde  data.  The  initial 
radius  was  2000  km,  and  it  was  decreased  by  400  km 
in  each  of  four  successive  passes  through  the  data.  The 
result  of  this  analysis  was  actually  worse  than  the  single 
scan  with  a  radius  of  750  km.  The  average  difference 
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24  hr.  ovBrog*  vs 
radii  of  influence 
Zonal  component 


Fig.  10.  (a)  Differences  between  the  24  h  average  radar  winds  and 
the  geostrophic  winds  for  various  tadii  of  influence  for  the  zonal 
wind  component. 


in  the  rms  error  was  slightly  greater  than  1  ms-1.  The 
successive  scans  method  produced  more  variability  in 
the  geostrophic  winds  calculated  from  the  analyzed 
geopotential  height  fields,  but  the  short-term  variations 
did  not  improve  the  agreement  with  the  radar  or  ra- 
winsonde  winds. 


24  hr  average  vs 
radii  of  influence 
Meridional  component 


Fig.  10.  (b)  As  in  Fig.  10a  but  for  the  meridional  component. 


8.  Correlation  between  winds  and  geopotential  heights 

The  correlation  coefficients  between  the  zonal  winds 
and  geopotential  heights  are  plotted  in  Fig.  11.  Mul¬ 
tivariate  objective  analysis  schemes  such  as  those  of 
Schlatter  (1975),  Rutherford  (1972)  and  Bergman 
(1979)  carry  out  the  objective  analysis  procedure  based 
on  both  height  and  wind  data.  The  wind  information 
supplements  the  observed  height  field  data  when  the 
two  quantities  are  related  through  the  geostrophic  re¬ 
lation.  The  procedure,  in  effect,  assumes  that  the  cor¬ 
relation  between  the  height  fields  and  the  winds  is 
unity  at  zero  spatial  separation.  Fig.  1 1  shows  the  cor¬ 
relation  coefficient  for  a  24  h  average  wind  to  be  nearly 
independent  of  height  with  a  value  of  approxi¬ 
mately  0.75. 

Since  the  statistics  for  the  rms  differences  are  clearly 
very  different  for  the  first  and  second  half  of  the  data 
sets,  the  correlation  coefficients  were  also  calculated 
separately  for  each  half.  The  results  are  shown  in  Table 
2.  The  correlations  are  nearly  identical  in  both  halves, 
unlike  the  situation  for  the  rms  differences.  When  the 
averaging  interval  is  increased  from  1  to  24  h,  the 
correlations  increase  for  both  the  radar/ rawinsonde 
comparison  and  the  radar/geostrophic  wind  compar¬ 
ison.  The  same  is  not  true  for  the  rawinsonde/geo- 
strophic  wind  comparison.  The  correlation  coefficients 
are  virtually  identical  for  the  l  h  average,  i.e.,  single 
profile,  and  24  h  average,  i.e.,  three-profile  compari¬ 
sons.  Again,  the  filtering  applied  to  the  height  variations 
in  the  balloon  wind  measurements  is  a  possible  ex¬ 
planation  for  this  result. 

9.  Comparison  with  earlier  results 

Most  previous  estimates  of  the  synoptic  error  of  the 
rawinsonde  wind  measurements  have  involved  com- 


Corrffotfons  for  24  tt  av«rogn 
Zonal  component 
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Fig.  11.  Cross  correlation  between  the  24  h  average  radar  and 
geostrophic  winds,  between  the  radar  and  rawinsonde  winds,  and 
between  averages  of  three  successive  rawinsonde  measurements  and 
the  geostrophic  winds.  The  correlations  are  for  the  zonal  wind  com¬ 
ponent  Values  for  the  meridional  component  are  similar. 
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Table  2.  Cross  correlations. 


3/1/79- 

3/17/79- 

3/16/79 

4/1/79 

u  u 

u 

V 

1-hour  averages 

Radar  vs.  Rawinsonde 

0.541  0.529 

0.522 

0.550 

Radar  vs.  Cressman  Analysis 
Rawinsonde  vs.  Cressman 

0.589  0.470 

0.545 

0.504 

Analysis 

0.661  0.734 

0.676 

0.829 

24-hour  averages 

Radar  vs.  Rawinsonde 

0.768  0.712 

0.702 

0.834 

Radar  vs.  Cressman  Analysis 
Rawinsonde  vs.  Cressman 

0.700  0.503 

0.682 

0.685 

Analysis 

0.691  0.731 

0.712 

0.826 

parisons  of  observed  values  and  the  output  grid  values 
from  numerical  models  or  objective  analysis  proce¬ 
dures  (e.g.,  Lyne  et  al..  1982;  Bergman,  1979;  Bengtsson 
et  al.,  1982).  Such  an  approach  was  necessitated  by 
the  lack  of  another  routine  measurement  technique 
with,  accuracy  comparable  to  the  rawinsonde  mea¬ 
surement  as  shown  by  Bengtsson  et  al.  (1982)  who 
present  a  summary  of  the  estimated  accuracies  of  the 
various  wind  measurements  used  operationally.  Also, 
the  comparison  between  observations  and  objective 
analysis  or  model  output  is  of  direct  interest  since  it 
provides  a  direct  estimate  of  the  magnitude  of  the 
unbalanced  wind  component  in  the  observations. 

Lyne  et  al.(\ 982)  found  the  difference  to  be  between 
3  m  s'1  at  low  levels  and  4  ms*1  at  high  levels  near 
100  mb.  Their  results  are  roughly  consistent  with, 
though  smaller  than,  the  range  of  2-6  m  s'1  given  by 
Bengtsson  et  al.  (1982).  Both  sets  of  differences  are 
much  smaller  than  the  best  case  estimate  of  7.98  m 
s"1  at  500  mb  calculated  by  Bergman  ( 1 979).  However, 
the  latter  was  based  on  assumed  values  for  the  mag¬ 
nitude  of  the  errors  that  are  needed  for  the  optimum 
interpolation  objective  analysis  scheme.  The  differences 
found  in  this  study  by  direct  comparison  of  the  balloon 
and  radar  measurements  were  3. 19  m  s'1  for  the  zonal 
component  and  4.03  m  s'1  for  the  meridional  com¬ 
ponent,  in  good  agreement  with  the  results  of  Lyne  et 
al.  (1982). 

Although  the  differences  between  the  two  wind 
measurements  are  likely  due  to  small-scale  spatial  and 
temporal  variability,  the  fact  that  the  difference  between 
the  measured  winds  and  the  geostrophic  winds  remains 
nearly  constant  as  the  averaging  interval  is  increased 
indicates  that  there  is  a  contribution  to  the  variance 
from  an  ageostrophic  wind  component.  The  rms  value 
of  the  ageostrophic  wind  component  based  on  the  dif¬ 
ferences  given  in  Table  1  would  be  approximately  1.5- 
2.0  ms'1.  This  study  has  dealt  with  only  one  site  and 
the  result  may  be  linked  to  the  location  and  local 
orography.  It  may  be  possible,  however,  to  sort  out 
such  effects  in  the  future  by  using  the  Poker  Flat  data 
sets  that  include  both  the  horizontal  and  vertical  wind 
components. 


Table  1  shows  that  there  is  a  clear  difference  between 
the  first  half  of  the  period  studied  and  the  second  half. 
In  the  first  half,  the  two  wind  measurements  are  very 
similar  and  the  estimate  of  the  ageostrophic  wind  com¬ 
ponent  is  the  same  based  on  both  techniques.  In  the 
second  half,  the  differences  between  the  radar  and  bal¬ 
loon  wind  measurements  are  much  larger  than  in  the 
first  half  and  the  difference  between  the  measured  and 
geostrophic  winds  is  about  the  same  as  the  difference 
between  the  two  independent  wind  measurements.  The 
data  processing  for  the  radar  data  did  not  change  in 
any  way  during  the  month  and  none  of  the  radar  com¬ 
ponents  such  as  transmitters  or  receivers  were  modified 
or  failed  during  this  period.  Thus,  there  is  no  reason 
to  expect  a  deterioration  in  the  quality  of  the  radar 
data.  Also,  the  comparisons  involving  the  balloon 
measurements  are  as  poor  at  this  time  as  those  in¬ 
volving  the  radar.  I  have  examined  the  500  mb  charts 
for  the  entire  period,  but  there  were  no  fundamental 
differences  in  the  overall  flow  pattern  during  the  first 
and  second  half  of  the  month.  One  difference  that  is 
evident  in  Fig.  1  is  the  amount  of  short-term  variability 
in  the  radar  data  after  the  time  labeled  600  h.  Aside 
from  the  episode  of  very  large  and  rapid  variations  at 
720  h,  there  is  also  consistently  more  variability  overall. 
Table  3  shows  the  variance  of  the  radar  winds  at  each 
height  for  the  first  and  second  half  of  the  period  and 
supports  this  conclusion.  It  is  clear  that  the  amount 
of  power  in  small  time-scale  variability  increases  in 
the  second  half  of  the  period,  although  the  cause  is 
not  clear  and  will  have  to  be  investigated  further. 

10.  Conclusion 

This  study  has  shown  that  a  VHF  radar  operated 
continuously  in  a  mode  such  as  the  Poker  Flat  MST 
radar  can  provide  synoptically  meaningful  data.  Com¬ 
parison  of  the  radar  wind  measurements  and  radio¬ 
sonde  pressure  and  wind  data  has  made  it  possible  to 
assess  the  representativeness  of  the  radar  data  directly. 
Differences  between  the  radar  and  balloon  wind  mea¬ 
surements  were  explained  by  the  wind  variability  over 
the  distance  separating  the  two  measurements  sites. 
As  the  averaging  interval  was  increased,  the  differences 
decreased. 

If  the  radar  data  are  used  in  the  future  for  synoptic 
applications,  the  data  will  either  have  to  be  averaged 


Table  3.  Variances  of  radar  winds. 
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temporally  over  some  period  appropriate  to  the  nu¬ 
merical  model  that  is  used  or  a  vertical  variation  filter 
such  as  that  applied  to  the  rawinsonde  data  will  have 
to  be  used.  The  results  of  this  study  will  help  in  de¬ 
signing  criteria  for  the  former  case.  The  latter  case 
could  not  be  investigated  using  Poker  Rat  data  that 
are  presently  available  since  the  height  resolution  is 
too  coarse  to  apply  the  empirical  rawinsonde  filters. 

The  data  set  studied  here  was  the  first  long  time 
series  of  horizontal  winds  produced  by  the  Poker  Flat 
radar.  Beginning  in  the  early  part  of  1981,  the  radar 
continuously  measured  all  three  wind  components  si¬ 
multaneously.  The  longer  data  set  should  be  used  to 
check  the  preliminary  results  presented  here,  including 
the  possibility  of  changes  in  the  short-term  variability. 
Finally,  the  same  types  of  comparisons  should  be  car¬ 
ried  out  with  some  of  the  better  objective  analysis 
schemes  that  are  now  available.  In  particular,  the  nor¬ 
mal  mode  method  has  many  advantages  in  that  it  not 
only  produces  analyzed  fields  that  are  exactly  consistent 
with  a  given  numerical  model,  but  it  can  also  be  used 
to  determine  the  contributions  from  the  fast  (gravity 
wave)  and  slow  (inertio-gravity  and  Rossby  wave) 
components  of  the  geopotential  and  velocity  fields  sep¬ 
arately.  Also,  the  normal  mode  analysis  does  not  im¬ 
pose  the  constraint  of  geostrophy. 

The  radar  measurements  have  been  shown  to  pro¬ 
duce  data  of  a  quality  that  is  at  least  comparable  to 
the  earlier  estimates  of  the  representativeness  of  the 
rawinsonde  winds  when  the  radar  data  are  averaged 
for  an  hour  or  more.  Of  course,  that  is  not  the  limit 
of  the  usefulness  of  the  radar  data  since  it  can  routinely 
provide  1  h  or  even  higher  time  resolution  data.  Such 
a  capability  will  become  increasingly  useful  as  the 
smallest  scales  resolved  by  numerical  models  decreases 
further. 
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1.  INTRODUCTION 

The  SOUSY-VHF-Radar  is  a  pulsed  coherent 
radar  operating  at  53.5  MHz  and  located  near 
Sad  Lauterberg,  West  Germany.  Since  1977,  the 
facility,  operated  by  the  Max-Planck-Institut 
fur  Aeronomie,  has  been  used  to  make  a  series 
of  frontal  passage  observations  in  the  spring 
and  fall.  Experiments  in  winter  have  been  dif¬ 
ficult  because  part  of  the  transmitting  and 
receiving  array  is  usually  covered  by  snow 
during  that  part  of  the  year.  Wavelengths 
around  6  m  are  known  to  be  sensitive  to  the 
vertical  temperature  structure  of  the  atmosphere 
(Green  and  Gage,  1980;  Rastogi  and  RSttger, 
1982).  Thus,  it  has  been  possible  to  use  radars 
operating  at  frequencies  near  50  MHz  to  locate 
the  tropopause.  Comparisons  between  radar  data 
and  radiosonde  data  have  shown  that  there  is  a 
large  gradient  in  the  radar  reflectivity  at  the 
height  where  the  radiosonde  tropopause  occurs. 

An  experiment  carried  out  by  Rottger  (1979) 
on  March  15-16,  1977,  showed  that  the  radar's 
sensitivity  to  the  vertical  temperature  struc¬ 
ture  could  also  be  used  to  locate  the  position 
of  fronts.  The  SOUSY-VHF-Radar  consists  of  a 
transmitting  array,  also  used  for  receiving  in 
some  configurations,  that  can  be  scanned  in  the 
off-vertical  direction  but  not  at  sufficiently 
low  elevation  angles  to  3tudy  the  horizontal 
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extent  of  structures  as  extended  as  fronts. 

Gage  and  Balsley  (1978),  Balsley  and  Gage  (1980), 
RSttger  (1980),  and  Larsen  and  RBttger  (1982) 
have  reviewed  UHF  and  VHF  Doppler  radar  techni¬ 
ques  and  applications  to  atmospheric  research. 

In  the  experiments  described  here ,  the  radar  was 
operated  in  the  spaced  antenna  mode.  Rottger  and 
Vincent  (1978)  and  Vincent  and  RSttger  (1980) 
have  described  the  method  and  its  advantages. 

The  transmitting  array  consists  of  196  Yagi 
antennas,  and  the  receivers  are  three  separate 
arrays  of  32  Yagls  each.  The  effective  antenna 
aperture  was  2500  m2,  the  applied  average  trans¬ 
mitter  power  was  typically  20  KW,  and  the  height 
resolution  was  150  m.  Vertical  profiles  of  the 
reflectivity  were  obtained  in  each  of  the  three 
receiver  arrays,  and  it  was  found  that  besides 
the  enhancement  of  the  signal  strength  associa¬ 
ted  with  the  tropopause  region,  there  was  also 
a  secondary  band  of  enhanced  reflectivities 
stretching  downward  from  the  upper  to  the  lower 
troposphere.  Comparisons  between  the  radar  data 
and  data  from  a  nearby  radiosonde  station  show 
that  the  band  Is  associated  with  the  cemperature 
gradients  In  a  passing  frontal  zone.  The  verti¬ 
cal  and  horizontal  velocities  were  also  measured 
during  the  experiments,  but  they  will  not  be  dis¬ 
cussed  here. 

The  results  of  analyzing  two  events  have  been 
presented  by  Rottger  (1979,  1981),  Rottger  and 
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Fig.  1.  Radar  reflectivities  as  a  function  of  height  and  time  measured  by  the  SOUSY-VHF-Radar. 
Observations  from  the  meteorological  observatories  at  Kassel,  60  km  southwest,  and  Hannover,  90  km 
north-northwest,  are  shown  on  the  left  and  right  sides  of  the  boxes,  respectively. 
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1981 


Fig.  2.  Radar  reflectivities  as  a  function  of 
height  and  time.  Contour  Interval  is  2  dB. 
Stippled  areas  show  higher  reflectivities. 


Schmidt  (1981) ,  and  Larsen  and  Rottger  (1983) . 
However,  in  this  paper  we  would  like  to  summa¬ 
rize  a  more  complete  data  set  consisting  of  a 
series  of  five  observations  of  frontal  struc¬ 
ture  made  with  the  SOUSY-VHF- Radar  in  March 
1977,  March  1981,  November  1981,  February  1982, 
and  April  1984.  The  extensive  data  set  shows 
results  essentially  in  agreement  with  Che  pre¬ 
liminary  results.  Comparison  of  cime/height 
cross  sections  of  reflectivity  measured  with  the 
radar  and  potential  refractivity  gradients  cal¬ 
culated  from  radiosonde  data  taken  at  a  nearby 
location  show  good  agreement.  Therefore,  we 
conclude  chat  the  radar  is  detecting  the  tem¬ 
perature  structure  associated  with  the  front  and 
that  Che  enhancement  in  reflectivities  is  not 
due  to  precipitation  or  other  very  localized 
processes.  Also,  we  have  found  that  the  radar 
can  be  used  to  locate  the  fronts  consistently, 
even  in  some  cases  when  the  fronts  are  rather 
weak. 


2.  March  15-16,  1977 

At  0000  UTC  on  March  15,  1977,  a  low  pres¬ 
sure  center  in  the  North  Atlantic  was  propaga¬ 
ting  eastward  toward  the  3rltish  Isles.  The 
associated  warm  front  and  the  trailing  cold 
front  extended  southward  from  the  center  of  the 
low.  The  warm  front  at  the  surface  had  traver¬ 
sed  France  and  West  Germany  by  0000  UTC  on 
March  17.  Reflectivities  measured  with  the 
SOUSY-VHF-Radar  from  1200  UTC  on  March  15  to 
0900  UTC  on  March  16  are  shown  in  Figure  1. 


Fig.  3.  Gradient  of  potential  refractive  index 
calculated  from  Hannover  radiosonde  data  for  the 
time  period  corresponding  to  the  reflectivities 
shown  in  Figure  2.  Stippled  areas  correspond  to 
larger  values. 

Cloud  cover  observations  from  two  nearby  meteor¬ 
ological  observatories  are  shown  below  the 
reflectivities. 

Of  paticular  interest  is  the  band  of  enhanced 
echoes  stretching  from  4  km  altitude  at  0000  UTC 
on  March  16  to  2  km  altitude  at  0600  UTC.  Extra¬ 
polating  the  slope  of  the  band  to  the  surface 
gave  a  time  for  the  surface  frontal  passage  in 
agreement  with  that  derived  from  the  weather 
charts.  Comparison  between  the  temperature  cross 
section  perpendicular  to  the  front  and  the  fea¬ 
tures  seen  in  the  radar  reflectivity  data  showed 
good  agreement . 


3.  March  6-7,  1981 

Figure  2  shows  the  reflectivities  observed 
with  the  radar  between  1200  UTC  on  March  6  and 
0900  UTC  on  March  7,  1981.  The  contour  interval 
is  2  dB,  and  the  stippling  indicates  regions  of 
higher  reflectivity.  The  signal  strength  gene¬ 
rally  decreases  with  height  in  the  troposphere 
but  increases  by  10-12  dB  over  a  vertical  dis¬ 
tance  of  0.5  km  or  so  ac  the  tropopause.  The 
tropopause  determinations  from  Hannover  radio¬ 
sonde  data  is  shown  by  the  heavy  bars  in  the 
figure.  The  bar  corresponding  to  0900  UTC  was 
actually  obtainded  at  1200  UTC,  the  scandard 
synoptic  time. 

The  feature  of  particular  interest  is  the 
band  of  enhanced  echo  strength  stretching  down¬ 
ward  from  the  tropopause  beginning  at  approxi¬ 
mately  1700  UTC.  The  feature  is  associated  with 
a  warm  front  chat  was  essentially  parallel  to  the 
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Fig.  4.  Reflectivities  measured  by  the  SOl'SY-VHF-Radar  during  the  period  from  Nov.  4  to  10,  1981. 


SNW-SSE  direction  and  propagated  eastward  across 
West  Germany  and  past  the  radar.  The  radar 
reflectivity  is  proportional  to  M2  with  M  given 
by 


M 


-77.6  x  10"6 
•{1+15|00 


P  r31n8, 

T  1  3z  ! 

( l  1  31 nq/ 3z.  , 
u  "  2  31n0/3z;/ 


(1) 


where  P  is  in  millibars,  T  is  absolute  tempera¬ 
ture,  3  is  potential  temperature,  and  q  is  the 
specific  humidity.  The  contours  of  M2  calcu¬ 
lated  from  the  Hannover  radiosonde  data  are 
shown  in  Figure  3  with  a  contour  interval  of 
4  dB.  Once  again,  the  stippled  areas  correspond 
to  higher  values.  The  agreement  between  the 
observations  and  the  calculated  values  is  quite 
good.  The  contours  of  M2  show  the  Increase  in 
the  height  of  the  tropopause  at  1200  UTC  on 
March  7.  Also,  the  enhancement  of  the  echoes 
between  6.8  and  8.3  km  after  0000  UTC  on  March  7 
is  apparent.  There  is  an  indication  of  the 
frontal  echo  band  in  the  values  of  the  potential 
refractivity  at  0000  UTC  as  shown  by  the  feature 
near  8.75  km  altitude.  However,  the  time  reso¬ 
lution  of  the  radiosonde  data  is  noc  sufficient 
to  show  the  details  of  the  frontal  zone  struc¬ 
ture. 


4.  November  4-10,  1981 

In  the  first  half  of  November  1981  observa¬ 
tions  were  made  over  a  period  of  more  than  a 
week.  The  reflectivities  for  the  period  from 
November  4  (Day  308)  to  November  10  (Day  314) 
are  shown  in  Figure  4.  On  Day  308,  the  reflec¬ 
tivities  decrease  with  altitude  in  the  tropo¬ 
sphere  and  then  begin  to  increase  Just  below  the 
tropopause.  A  cold  frontal  band  is  observed 
after  1600  UTC  and  stretches  upward  in  altitude 
with  time  as  would  be  expected.  Two  warm  frontal 
bands  are  evident  as  enhanced  reflectivity 
regions  moving  downward  with  time.  The  first  is 
observed  near  the  tropopause  at  0000  UTC  on 
Day  309.  The  second  is  first  seen  at  0600  UTC 
on  Day  310 . 

Potential  temperatures  were  calculated  from 
the  Hannover  radiosonde  data  and  have  been  plot¬ 
ted  and  contoured  on  a  scale  similar  to  that  used 
for  the  reflectivities.  The  results  are  shown  in 
Figure  5.  The  same  general  features  are  evident 
in  both  figures.  The  packing  of  the  potential 
temperature  contours  is  characteristic  of  the 
tropopause  and  stratosphere.  The  lower  boundary 
of  the  packing  is  located  at  the  same  height  as 
the  increase  in  the  radar  reflectivities  associa¬ 
ted  with  the  tropopause.  The  cold  front  and  the 


Julian  Day  1981 


Fig.  5.  Potential  temperatures  calculated  from  the  Hannover  radiosonde  data  for  the  period  correspond¬ 
ing  to  the  radar  reflectivities  shown  in  Figure  4. 


7  FEB  12 


00  8  FEB  12 


00  9  FEB  82  12  UTC 


tig.  6.  Same  as  Figure  2  but  for  Feb.  7-9,  1982. 

second  warm  front  show  up  clearly  in  the  poten-  5.  February  7-9,  1982 

tial  temperature  data,  but  the  first  warm  front¬ 
al  band  is  not  as  clearly  evident.  However,  The  reflectivities  measured  by  the  radar 

satellite  photos  for  this  time  show  two  distinct  during  the  period  from  February  7  to  9,  1982, 

cloud  bands.  are  shown  in  Figure  6.  Once  again  the  stippled 

During  the  time  form  Nov.  4  to  10,  the  radar  areas  represent  larger  reflectivities.  However, 

site  was  in  a  region  of  northerly  flow  on  the  this  time  the  intermediate  values  have  been  cross 

eastern  side  of  a  stationary  high  pressure  hatched,  as  well.  The  typical  features  already 

system  centered  over  the  British  Isles.  The  noted  in  the  previous  examples  are  present.  The 

fronts  that  traversed  the  radar  site  were  con-  echo  strength  decreases  with  altitude  in  the 

fined  to  the  upper  troposphere,  as  the  potential  troposphere,  but  there  is  an  enhancement  in  the 

temperature  cross  section  shows.  reflectivities  in  the  upper  troposphere  associa¬ 

ted  with  the  temperature  discontinuity  that  de¬ 
fines  the  location  of  the  tropopause.  The  heavy 
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Fig.  7.  Same  as  Figure  3  ouc  corresponding  to  the  data  shown  in  Figure  8. 


Fig.  8.  Same  as  in  Figure  2  but 

circles  show  the  tropopause  height  determined 
from  the  Hannover  radiosonde. 

Near  1200  UTC  on  February  8,  a  band  of 
enhanced  reflectivities  associated  with  a  warm 
front  stretches  downward  from  the  tropopause. 
Beginning  at  approximately  0200  UTC  on 
February  8,  a  band  of  intermediate  reflectivi¬ 
ties  stretches  upward  in  connection  with  a  cold 
front  associated  with  the  low  pressure  system. 

The  cold  fronc  for  chis  particular  event  was 
much  weaker  than  the  warm  front,  and  the  dif¬ 
ference  in  strength  accounts  for  the  difference 
in  che  magnitude  of  the  reflectivities. 

We  have  also  calculated  the  potential  re¬ 
fractive  index  gradient  M2  based  on  Hannover 
radiosonde  data.  The  vertical  time  section  is 
shown  in  Figure  7.  The  same  features  are  evi¬ 
dent  in  both  the  reflectivity  and  M2  cross 
sections.  More  detailed  analysis  of  the  event 
is  given  by  Larsen  and  Roctger  (1983). 


6.  April  9-13,  1984 

Observations  were  made  with  the  S0USY-VHF- 
Radar  over  a  five  day  period  in  April  1984. 

From  1200  UTC  on  April  9  until  1200  UTC  on 
April  11,  che  radar  was  located  in  a  region  of 
northerly  flow  on  the  eastern  side  of  a  sta¬ 
tionary  high  pressure  system.  During  this  time 
a  number  of  mesoscale  disturbances  with  hori¬ 
zontal  scales  of  50  to  100  km  developed  in  the 
region.  Most  did  not  pass  the  radar,  but  one 
such  system  was  observed  near  0000  UTC  on 
April  11.  The  radar  reflectivities  are  shown 
in  Figure  3. 

The  signal  strength  showed  an  enhancement 
in  the  upper  troposphere  during  che  passage  of 
the  mesoscale  system,  but  the  feature  is  not  as 
narrow  and  does  not  show  che  tilt  characteristic 
of  the  ocher  frontal  passage  observations  chat 
we  have  presented.  Lace  in  the  day  on  April  12 
a  surface  froncal  passage  took  place,  as  shown 
by  the  weather  charts  and  satellite  photographs. 
The  features  typical  of  a  frontal  passage  are 
clearly  present  in  che  reflectivity  data  begin¬ 
ning  at  about  1600  UTC  on  April  11.  A  band  of 
enhanced  echoes  begins  to  descend  from  che 
height  of  che  tropopause  and  screeches  downward 
toward  the  surface.  The  air  mass  following  the 
froncal  passage  then  has  a  higher  tropopause 
height  chan  the  air  mass  prior  to  the  frontal 
passage.  The  morphology  in  chis  event  is  very 
similar  to  the  morphology  in  che  other  four 
events  presented  here. 


7.  Conclusion 

Observations  of  five  separate  froncal  pas- 


for  the  period  April  9-13,  1984. 

sage  events  made  with  the  SOUSY-VHF-Radar  have 
shown  that  a  radar  operating  at  a  frequency  near 
50  MHz  can  be  used  to  detect  the  location  of 
fronts  on  a  routine  basis.  The  SOUSY-VHF-Radar 
does  not  operate  on  a  continuous  basis.  While 
the  observation  periods  were  chosen  on  the  basis 
that  it  would  be  likely  that  a  frontal  passage 
would  take  place  at  the  location  of  the  radar,  no 
attempt  was  made  to  choose  only  the  strongest  or 
most  developed  fronts.  In  fact,  the  front  obser¬ 
ved  in  April  1984  and  the  cold  front  observed  in 
February  1982  did  not  have  particularly  strong 
temperature  gradients. 

The  comparison  between  radar  and  radiosonde 
data  have  shown  good  agreement,  indicating  that 
the  features  seen  in  the  radar  reflectivity  daca 
are  characteristic  of  che  froncal  temperature 
structure  and  are  not  associated  with  precipi¬ 
tation  or  local  convection.  That  is  particular 
true  since  the  separation  between  the  radar  site 
and  the  nearest  radiosonde  station  is  approxi¬ 
mately  90  km. 
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VHF-Radar  Observations  o  f  Frontal  Zone 
Structure 

M.  F.  Larsen  (Dept,  of  Physics  and  Astronomy, 
Clerason  University,  Clemson,  SC  29631)  and 
J.  Rottger  (Arecibo  Observatory,  P.  0. 
Box  995,  Arecibo,  PR  00612) 

Since  March  1977,  the  SOUS Y- V HF-Ra d a r  has 
been  used  to  make  a  series  of  frontal  passage 
observations  in  the  spring  and  fall.  The 
radar  operates  at  a  wavelength  of  5.6  m  and, 
thus,  is  sensitive  to  the  temperature  struc¬ 
ture  of  the  atmosphere.  It  has  already  been 
shown,  for  instance,  that  a  VHF  radar  can  be 
used  to  determine  the  location  of  the  tropo- 
pause.  We  present  data  from  five  separate 
experiments  showing  that  the  reflectivity  -data 
can  also  be  used  to  determine  the  location  of 
tropopause  breaks  and  frontal  boundaries  be¬ 
tween  the  surface  and  the  tropopause.  In  each 
case,  we  compare  reflectivities  measured  with 
the  radar  and  the  potential  refractivity  gra¬ 
dient  calculated  from  radiosonde  data  from  a 
station  90  km  away.  There  is  good  agreement 
between  the  calculated  and  observed  values, 
although  the  radar  data  show  more  of  the  de¬ 
tailed  structure  due  to  the  much  higher  time 
resolution . 
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This  abstract  describes  some  of  the  later  work  supported  by  AFOSR-83-lOO.  The 
paper  has  not  been  delivered,  in  fact,  but  the  work  summarized  here  will  be  the 
basis  for  a  future  publication. 


Short  Abstract 

OBSERVATIONS  OF  MESOSCALE  WAVE  STRUCTURE  WITH  A  VHF 

DOPPLER  RADAR 

i  n 

M.  F.  Larsen  and  J.  Rottger 

^Dept.  of  Physics  and  Astronomy 
Clerason  University 
Clemson,  SC  29631 

^Arecibo  Observatory 
P.  0.  Box  995 
Arecibo,  PR  00612 


Abstract.  In  November  1981  and  April  1984,  the  SOUSY-VHF-Radar 
located  in  West  Germany  operated  for  periods  of  16  days  and  5 
days,  respectively,  with  the  goal  of  obtaining  data  during  one  or 
more  frontal  passages.  In  both  cases,  mesoscale  wave  structure 
was  observed  when  weak  pressure  and  temperature  gradients 
dominated  at  the  surface  and  aloft.  The  waves  appeared  as 
downward-sloping  bands  of  enhanced  reflectivities  in  the  altitude 
range  from  0.45  km,  the  lowest  radar  range  gate,  and  up  to  the 
tropopause.  Three  to  five  oscillations  were  evident.  In 
November  1981,  the  wave  period  in  the  earth-fixed  reference  frame 
was  approximately  9  hr,  while  in  April  1984,  it  was  closer  to 
5  hr.  The  vertical  velocities  measured  with  the  radar  wind 
profiler  showed  30  to  40  cm/s  oscillations  with  the  same  vertical 
and  temporal  structure  as  the  r-e  f  1  e  c  t  i  v  i  t  y  bands,  although  90° 
out  of  phase.  The  satellite  photographs  for  both  periods  showed 


banded  cloud  structure  above  the  radar. 
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